AD- A  12 1  220 

UNCLASSIFIED 


I 


ON  DERIVING  AND  SOLVING  THE  GENERALIZED  BIVARIATE 
LINEAR  LOCATION  PROBLEMS(U)  AIR  FORCE  INST  OF  TECH 
WRIGHT-PATTERSON  AFB  OH  J  W  COLEMAN  SEP  82 
AF I T/C I /NR/82-610  _  _ 


NL 


F/G  12/1 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BuPfAU  0*  STANDARDS- }%*  A 


SECURITY  CLASSIFICATION  OF  ThiS  PAGE  (*1>en  D*!m  fntnrd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3  REClPlF*  T’S  CATALOG  NUMBER 


t.  REPORT  NUMBER 

AFIT/C I/NR/82-61 D 


♦.  TITLE  (*nd  Submit)  S.  TYPE  OF  FJ  IPORT  A  PERIOD  COVERED 

On  Deriving  and  Solving  the  Generalized  Bivariate,  ftf^^/DISSERTATION 
Linear  Location  Problems  _ 

fit  6  performing  org.  report  number 


7.  AuTHORf  •) 


Joseph  William  Coleman 


0.  CONTRACT  OR  GRANT  NUMBER(*j 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

AFIT  STUDENT  AT:  Arizona  State  University 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  6  WORK  UNIT  NUMBERS 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

AFIT/NR 
WPAFB  OH  45433 


12.  REPORT  DATE 


September  1982 


II.  NUMBER  OF  PAGES 

102 


14.  MONITORING  AGENCY  NAME  S  ADDRESS^//  dilferent  From  Controlling  Office)  IS.  SECURITY  CL  ASS.  (ol  thi»  report) 

UNCLASS 


I  17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


APPROVED  FOR  PUBLIC  RELEASE:  IAW  AFR  190-17 


29  OCT  1982 


1 19.  KEY  WORDS  (Continue  on  reveree  tide  it  necaeeary  and  identify  by  block  number) 


I  20.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 


ATTACHED 


lM  e.  WOLAVER  I 

Dean  for  Research  and  » 

Professional  Development 
AFIT.  Wriqht-Patterson  AFB  lOH 


92  1  1  09  09  S 


DD  1  j an*73  1473  EDITION  OF  1  NOV  «S  IS  OBSOLETE 


UNCLASS 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fRTien  D»lm  Enffd) 


ABSTRACT 


The  generalized,  bivariate,  linear  location  problem 
concerns  the  locating  of  a  linear  facility,  X2*=Cq+C^x^, 
in  a  two-dimensional  Euclidean  space  such  that  the  p-norm 
distance  taken  to  the  q  power  is  minimized  of  serving  n 
existing  fixed  facilities  whose  location  in  the  two- 
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dimensional  Euclidean  space  xs  given  by 
i=l,...,n.  The  solution  of  the  generalized,  bivariate, 
linear  location  problem  consists  of  two  subproblems.  The 
first  subproblem  involves  the  determination  of  the  point 
on  any  linear  facility  that  minimizes  the  p-norm  distance 
to  an  individual  existing  facility.  The  second  subproblem 
consists  of  determining  the  optimal  linear  facility  that 
minimizes  the  sum  of  the  q  multiples  of  the  p-norm  dis¬ 
tance  from  all  the  existing  facilities  to  the  point  on  the 
linear  facility  determined  by  the  previous  step. 

The  lack  of  convexity  of  the  generalized,  bivari¬ 
ate  linear  location  problem  prohibits  a  universal  solution 
technique  for  all  combinations  of  possible  values  for  p 
and  q.  For  certain  combinations  of  p’s  and  q's,  however, 
an  exact  solution  can  be  determined.^  For  example,  the 


r  ■  if W«  - 


case  where  the  limit  is  taken  as  p  approaches  one  from 
the  negative  direction  and  q  equals  two  reduces  to  simple 


regression  of  on  x^.  If  the  limit  is  taken  as  p 
approaches  one  from  the  positive  direction  and  q  equals 
two,  the  generalized,  bivariate,  linear  location  problem 
becomes  simple  regression  of  on  x^.  In  general,  most 
combination  of  p's  and  q's  greater  than  one  can  only  be 
solved  by  a  heuristic  approximation  procedure. 
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CHAPTER  I 


INTRODUCTION 

The  generalized,  bivariate,  linear  location  problem 
concerns  the  placement  of  a  linear  facility  of  the  form 
x_  =  C  +  C.x,  such  that  certain  criteria  dictated  by  the 
decision  maker  are  minimized.  The  linear  facility  to  be 
located  might  be  considered,  for  example,  from  the  classi¬ 
cal  location-theory  perspective  to  be  a  pipeline  connect¬ 
ing  a  number  of  factories,  a  trunk  line  within  a  plant 
connecting  work  stations,  or  an  interstate  highway  con¬ 
necting  various  cities.  A  more  generic  interpretation  of 
the  linear  facility  might  involve  the  estimating  (or  fix¬ 
ing)  of  a  relationship  between  output  of  a  machine  and 
power  consumed  by  that  machine,  sales  of  a  product  and 
the  advertising  expenditures  on  that  product,  or  (even  in 
the  classical  sense)  the  size  of  the  seeds  of  the  daughter 
sweet  pea  plaints  with  respect  to  the  size  of  the  seeds  of 
the  mother  sweet  pea  plants  (David,  1978) . 

The  critical  aspect  of  the  generalized,  bivariate, 
linear  location  problem  concerns  the  criteria  selected 
by  the  decision  maker  to  be  minimized.  If  a  linear 
facility  cannot  pass  exactly  through  each  data  point,  the 


2 


person  performing  the  location  decision  must  determine 
which  distance  should  be  minimized.  Should  the  objective 
function  of  the  minimization  problem  consider  just  the 
distance  with  respect  to  one  of  the  variables,  or  should 
the  distance  to  be  minimized  be  some  composite  distance 
with  respect  to  both  variables?  In  addition,  the  decision 
maker  must  determine  the  relative  importance  of  each  of 
the  existing  facilities  with  respect  to  the  distance  being 
minimized.  For  example,  should  the  importance  of  an  exist¬ 
ing  facility  in  the  minimization  process  be  linearly  pro¬ 
portional  to  the  distance  to  the  linear  facility,  or 
should  outliers  have  more  or  less  importance? 

Early  contributors  to  point- location  literature 
favored  the  Euclidean  norm  of  modeling  distances  in  the 
minimization  process.  Later,  the  rectangular  norm  was 
used  for  approximating  distances  when  movement  was 
restricted  to  a  network  which  was  basically  a  rectangular 
grid  (Morris,  1981) .  Both  norms  are  special  cases  of  the 
p-norm  distance  function, 

Vx)  ■  1  -  ait!P’1/P'  P  i 

where  x  and  a^  are  points  in  N-dimensional  space.  Love 
and  Morris  (1972) ,  considering  the  point-location  problem, 
found  that  p-norm  distances  raised  to  the  power  of  p  =  1.69 
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and  1.78  were  better  models  for  interstate  road  distances 
them  were  the  classical  rectilinear ,  p  =  1,  and  Euclidean, 
p  =  2,  distances.  Cooper  (1968),  also  working  with  the 
point-location  model,  observed  that  the  weighting  of  vari¬ 
ous  data  points,  regardless  of  the  p-norm  distance  used, 
should  not  be  necessarily  directly  proportional  to  the 
distance  from  the  point  to  the  linear  facility.  Depending 
upon  the  situation,  the  p-norm  distance  taken  to  powers 
other  than  one  may  more  accurately  model  the  situation  in 
terms  of  costs.  Considering  the  location  (or  estimation) 
of  linear  models,  the  least  squares  procedure  used  in 
regression  minimizes  the  sum  of  the  squared  vertical 
distances  between  the  data  points  and  the  equation 
x2  *  CQ  +  C-jX^  (Neter  &  Wasserman,  1974);  and  the  least 
sum  procedure  minimizes  the  sum  of  the  absolute  values  of 
the  same  distances. 

The  generalized,  bivariate,  linear  location  problem 
allows  the  decision  makers  first  to  select  the  p-norm 
distance  taken  to  any  power  q  that  will  best  satisfy  the 
needs  of  their  analysis.  The  linear  facility  can  then  be 
located  using  the  data  with  respect  to  the  predetermined 
values  of  p  and  q. 

Statement  of  the  Problem 

The  generalized,  bivariate,  linear  location  problem 
can  be  stated  in  the  following  manner:  Given  n  existing 


fixed  facilities  i  whose  locations  in  a  two-dimensional 
Euclidean  space  are  given  by  (a^,  a ^2)  with  i  =  1,  ...»  n; 
locate  a  linear  facility,  x2  -  CQ  +  cixi'  the  two“ 

dimensional  Euclidean  space  such  that  the  total  cost  of 
serving  the  n  fixed  facilities  i  is  minimized.  Assume 
that  cost  is  proportional  to  distance  so  that  the  cost  of 
serving  facility  i  is  equal  to  where  WL  is  a  weight¬ 

ing  constant  which  transforms  distance  into  cost  and 

pi 

is  the  p-norm  distance  taken  to  the  q  power  from  facility 
i  to  the  linear  facility  which  is  to  be  located. 

More  specifically,  locating  a  linear  facility  con¬ 
sists  of  two  subproblems:  (1)  determine  the  point 
(x^,  x^2)  on  the  linear  facility  closest  to  the  existing 
facility  i,  and  (2)  determine  estimates  for  the  linear 
parameters,  CQ  and  C^. 

Formally,  the  modified  distance  function,  £j^(Co,  C^) , 
is  given  by 


l*. 

Pi 


<Co' 


Cl> 


2 

[  I 
t-1 


‘it 


Pjq/p 


subject  to 


where  x^t  is  the  tth  coordinate  of  a  point  (x^,  xi2^  on 
the  linear  facility,  x2  =  CQ  +  C^x^,  that  minimizes  the 
p-norm  distance  raised  to  the  q  power  from  the  ith  exist¬ 
ing  facility  to  the  linear  facility  to  be  located.  The 
preceding  definition  is  similar  to  the  work  of  Morris 
(1981)  with  the  addition  of  the  linear  constraint. 

Using  the  preceding  definition,  the  formal  statement 
for  the  linear  Weber  problem  would  be 


Minimize  Cl^lC  , 
p  o 


Cl>  * 


n 

Z 

i=l 


V?i(Co' 


ci> 


where 

CQ,  are  the  intercept  and  slope  parameters 

of  the  linear  facility  to  be  determined, 

is  the  weighting  constant  which  transforms 

distance  into  costs,  and 

(C  ,  C. )  is  the  minimum  p-norm  distance  from 
pi  o  l 

the  point  (ai;L»  ai2)  to  the  point 
(xil,  xi2)  on  the  line,  x2  =  CQ  +  cix]/ 
taken  to  the  power  q. 

Review  of  Related  Literature 
The  subject  of  generalized,  bivariate,  linear  location 


is  a  child  with  many  fathers.  If  the  problem  of  locating 
a  "linear"  facility  is  considered  an  extension  of  the 
point  location  area,  then  the  roots  of  the  linear  location 
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problem  may  be  traced  to  Fermat  in  the  early  17th  century 
(Wesolowsky,  1973) .  If  the  linear  location  problem  is 
considered  a  generalization  of  the  least  square  estimation 
procedure/  then  the  heritage  of  the  linear  location  problem 
lies  in  the  earliest  works  of  least  square  estimation  done 
by  Gauss  in  the  late  18th  century  (Eisenhart,  1978) .  To 
develop  a  proper  perspective  regarding  the  linear  loca¬ 
tion  problem,  therefore,  consideration  must  be  given  to 
both  location  theory  and  least  sum  regression  theory. 

Location  Theory 

Fermat  first  stated  the  point  location  problem  in 
the  following  manner:  "Given  three  points  in  the  plane, 
find  a  fourth  point  such  that  the  sum  of  its  distances  to 
the  three  given  points  is  a  minimum"  (Wesolowsky,  1973, 
p.  96).  Torricelli  found  the  point  that  minimized  the 
sum  of  its  distances  to  the  points  on  a  triangle  to  be 
"at  the  intersection  of  the  circles  which  circumscribe  the 
equilateral  triangles  constructed  on  the  sides  of  and  out¬ 
side  the  triangle  formed  by  the  three  given  points" 
(Verdini,  1976,  p.  1).  From  the  work  of  Cavalieri, 

Simpson,  and  Heinen,  Jacob  Steiner  posed  and  solved  the 
Steiner  problem. 

If  all  angles  of  the  triangle  formed  by  the  given 

points  are  less  than  120°,  then  the  minimizing  point 

is  found  at  the  point  where  three  lines,  one  from 


- 


each  given  point,  intersect  at  120°.  If  one  of  the 
angles  of  the  given  triangle  is  greater  than  or  equal 
to  120°,  the  vertex  of  this  angle  is  the  minimizing 
point.  (Verdini,  1976,  p.  1) 

The  distances  referred  to  in  these  solutions  were  all 

Euclidean  distances  raised  to  the  first  power  (p  =  2, 

q  =  1) ,  and  the  points  were  all  of  equal  weight.  Simpson, 

in  1750,  was  the  first  to  generalize  Fermat's  problem  by 

investigating  the  sum  of  weighted  distances  (Wesolowsky, 

1973)  . 

In  1909  Alfred  Weber  published  his  pioneering  work 
in  special  economics,  liber  den  Standort  der  Industrien. 

He  posed  the  problem  of  placing  a  factory  that  produced 
one  output  manufactured  using  two  raw  inputs  from  separate 
localized  areas  and  subsequently  sold  at  one  distant 
market.  Weber  assigned  the  three  elements  i  (i.e.,  one 
product  and  two  raw  inputs)  weights  u^  depending  upon  the 
weight  of  the  element  divided  by  the  weight  of  the  product 
With  this  system  the  weight  of  the  product  was  defined  as 
one.  Weber  further  defined  the  transportation  rates  for 
each  product  as  r^.  Thus,  "Weber  termed  'ideal  weights,' 
or  weights  of  elements  adjusted  for  differential  transpor¬ 
tation  rates"  (Kuhn  &  Kuenne,  1962,  p.  22)  to  be 
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The  total  transportation  cost,  C£p(x^,  x2>  to  be  minimized 
is  the  following: 


x2)  ~  i^1WiZ2i 


where  £2^  is  the  Euclidean  distance  from  each  facility 

that  supplies  the  raw  inputs  or  receives  the  final  product 

(ail»  ^i2)  to  the  plant  to  be  located  (x^  x2>  .  The 

Euclidean  distance  is  mathematically  defined  as 

2  2  1/2 

i-2i  *  Hx^  -  aii^  +  ^x2  =  ai2*  (Kuhn  &  Kuenne, 

1962) . 

The  earliest  solutions  to  the  Weberian  problem  were 
either  geometrical  or  physical  in  nature.  Such  models  as 
the  weighted  triangle  (Dean,  1938),  Launhart's  pole 
principle  (1885),  isodapones  (Palander,  1935),  and  the 
Varignon's  frame  (Isard,  1956)  provided  exact  solutions, 
although  they  were  quite  cumbersome  and  limited  in  the  size 
of  problem  that  they  could  solve  (Kuhn  &  Kuenne,  1962). 

Weiszfeld  was  the  first  to  discover  an  iterative 
technique  to  solve  the  original  Weberian  problem.  The 
technique  of  iteratively  solving  two  simultaneous  equa¬ 
tions,  first  published  in  1937,  was  subsequently  redis¬ 
covered  in  the  fifties  and  sixties  by  various  authors 
(Cooper,  1963;  Kuhn  &  Kuenne,  1962;  Michle,  1958).  To 
minimize  the  cost  function,  the  first  derivatives  are  taken 


with  respect  to  x^  and  x2,  and  the  results  are  set  equal 
to  zero.  This  process  yields 


3C&2 (X) 

n 

—  v 

VX1  - 

3x^ 

w 

i=l 

»2iUl, 

X2> 

3C&2  (X) 

n 

—  y 

Vx2  ' 

ail) 

3x2 

i=l 

*2i(xl' 

X2> 

Rearranging  terms  produces  the  following  equations: 

"  Wiail 

_  i=l  *2i(xl'  x2] 

X1  "  n  Wi 

i=l  *2] 

?  Vi2 

-  i=l  £2i<xl*  x2} 

2  ”  n  VT 

i*l  42i*xl'  x2* 


An  iterative  procedure  is  required  since  x.^  =  f(x^,  x2> 
and  x2  =  f(x1#  x2).  The  weighted  mean,  the  center  of 
gravity,  was  a  convenient  starting  point  for  the  iteration. 
The  preceding  function  was  found  to  be  convex  since  it  was 
the  sum  of  convex  functions,  and  the  iterative  procedure 
was  assumed  to  converge  on  the  global  minimum  if  no  itera¬ 
tion  resulted  in  x^  and  x2  being  equal  to  the  coordinates 
of  one  of  the  existing  facilities. 
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Cooper  (1968)  generalized  the  Weber  problem  to  mini¬ 
mize  the  Euclidean  distance  raised  to  a  power  q  where  q 
was  greater  than  zero.  Love  and  Morris  (1972)  extended 
the  Weber  problem  to  include  p-norm  distances,  but  Morris 
(1981)  combined  these  two  formulations  into  the  generalized 
Weber  problem  which  can  be  stated  in  the  following  manner 
for  Euclidean  two  space : 


n 

Minimize  ,  x_)  »  Z  W.Jl^.fx,  ,  x_) 

Pi  2  i_i  x  pi  1  2 


i=l 


where 


*pi(Xl'  X2>  = 


2 

l  Z  I 
t=l 


it 


'  ait 


?] q/p . 


Inspired  by  the  work  of  Love  and  Morris  (1972),  Wesolowsky 
and  Love  (1972),  Eyster,  White,  and  Wierwille  (1973), 
and  Love  (1969),  Verdini  (1976)  developed  and  proved  con¬ 
vergence  of  a  generalized  hyperbolic  approximating  func¬ 
tion  that  eliminated  the  difficulties  in  differentiating 
the  generalized  Weber  problem  for  q  -  1  and  p  =  1  or  2. 

The  approximating  function,  X£),  took  the  follow¬ 

ing  form: 


Lpi<x 


1' 


X2} 


=  { 


2 

I 

t=l 


l(Xit 


‘  ait) 


2  +  £]p/2}q/p 
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where  e  was  a  strictly  positive  smoothing  constant. 

Morris  and  Verdini  (1979)  extended  the  application  of  the 
hyperbolic  approximating  function  by  proving  convergence 
and  convexity  for  q  =  1  and  p  >_  1,  and  Morris  (1981) 
further  generalized  the  convexity  and  convergence  proofs 
for  q  >.  1  and  p  >  0.  As  Love  and  Morris  (1978)  proved  for 
p  =  2  and  Morris  (1981)  proved  for  p  >  0,  however,  the 
Weber  problem  is  not  convex  or  concave  for  q  <  1.  Local 
minima  exist  at  each  existing  facility. 

Location  theory  was  extended  from  the  point  perspec¬ 
tive  to  the  linear  perspective  by  MacKinnon  and  Barber 
(1972)  when  they  developed  an  algorithm  that  positioned 
a  linear  facility  that  minimized  the  Euclidean  distance 
(p  =  2,  q  =  1) .  They  used  a  heuristic  method  to  search 
for  the  optimal  slope  of  the  linear  facility  between 
bounds  of  the  slope  of  the  regression  lines  of  x2  on  x^^ 
and  x^  on  x2. 

Wesolowsky  (1975)  determined  an  exact,  although  not 
necessarily  an  unique,  solution  to  the  linear  location 
problem  for  Euclidean  distance  (p  =  2,  q  =  1)  by  rotating 
the  original  coordinate  axes  until  the  optimal  slope  of 
the  linear  facility  was  parallel  to  the  transformed  axes. 
The  points  (x^,  xi2^  where  i  =  If  ...»  n  were  transformed 
into  polar  coordinates  (x^lp,  xi20*  the  f°llowin9 


manner : 
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Xi29  = 


x 


iie 


x^2  cos  9  +  x^  sin  9 , 
-x^2  sin  9  +  cos  9  . 


Thus,  with  the  knowledge  that  the  line  that  minimizes  the 
p  =  2,  q  =  1  distance  is  the  median  line  (Rice  &  White, 
1964) ,  Wesolowsky  found  the  value  of  0  for  the  line 
xi20  =  Z Q  in  the  transformed  axes  that  minimized  the 
following  equation: 


Cl^{Q)  = 


n 
Z  W 
i=l 


i9 


Z9  " 


xiie 


Once  the  value  of  0  was  found,  was  equal  to  tan  6;  and 

C  was  the  median  of  x.OQ. 
o  i 

Morris  and  Norback  (1980)  simplified  Wesolowsky' s 
procedure  by  recognizing  that  the  linear  facility  for  q  =  1 
was  the  median  line  and  that  at  least  one  solution  would 
pass  through  at  least  two  existing  facilities  i.  Morris 
and  Norback  also  extended  the  applicability  of  their 
procedure  to  0  <  q  <_  1  distances.  The  authors  pointed 
out,  hoyever,  that  the  convexity  condition  for  q  -  1  was 
not  present  for  1  <  q  <  00 ,  and  that  an  alternate  heuristic 
technique  must  be  used  for  solving  problems  with  values  of 
q  >  1. 
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Least  Sum  Regression  Theory 

The  least  square  estimation  procedure,  which  is  the 
basis  for  most  linear  regression  applications,  was  first 
published  by  Gauss  in  1794.  His  treatment  of  the  location 
of  a  function  in  the  form  y  =  a  +  BX  which  minimized  the 
squared  error  from  a  set  of  points  n  where  n  >  2  was 
inspired  by  John  Neinrick  Lambert's  approach  to  the  same 
problem  by  using  a  moving  average  (Eisenhart,  1978) . 
Francis  Galton  indirectly  coined  the  term  "regression" 
in  his  1885  publication,  Natural  Inheritance,  when  he 
studied  sweet  pea  plants.  Galton  noticed  that  the  size 
of  the  seeds  of  daughter  plants  seemed  to  be  related  to 
the  size  of  the  seeds  of  the  mother  plant,  but  that  they 
also  appeared  to  "revert"  to  the  mean.  The  term  "revert" 
was  soon  replaced  by  "regress"  (David,  1978). 

A  relatively  short  time  after  Gauss,  in  the  1820s 
Fourier  first  proposed  a  method  of  minimizing  the  sum  of 
absolute  deviations  instead  of  minimizing  the  sum  of  the 
squared  deviations  as  used  in  regression.  Fourier  recom¬ 
mended  am  iterative  technique  that  he  developed  which  was 
similar  to  the  simplex  method  (Fisher,  1961) .  Later  in 
the  century,  Edgeworth,  receiving  inspiration  from  the 
works  of  Laplace,  recognized  that  the  least  squares  pro¬ 
cedure  was  highly  susceptible  to  outliers  and  was,  there¬ 
fore,  inferior  to  least  sum  estimators  (Barrodale,  1968). 
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This  observation  was  later  supported  by  the  studies  of 
Barrodale  (1968)  and  Forsythe  (1972).  The  algorithm  pro¬ 
posed  by  Edgeworth  (Rhodes,  1930) ,  as  well  as  more  recent 
algorithms  developed  by  Rhodes  (1930) ,  Singleton  (1940) , 
and  Karst  (1958) ,  were  either  limited  to  two  dimensions 
or  became  unwielding  as  the  number  of  dimensions  increased 
(Fisher,  1961) . 

Each  of  the  preceding  authors  recognized  that  mini¬ 
mizing  the  sum  of  the  absolute  deviations  did  not  neces¬ 
sarily  produce  a  unique  solution.  Rice  and  White  (1964) 
pointed  out  that  the  y-intercept  (Cq)  was  the  median  of 
ai2  -  Ciail;  and,  therefore,  the  best  least  stun  estimator 
would  be  the  median  line.  Sposito  and  Smith  (1976), 
expanding  on  the  work  of  Appa  and  Smith  (1973),  provided 
the  following  conditions  for  the  location  of  a  hyperplane 
of  R  dimensions  from  n  observations  which  minimized  the 
sum  of  absolute  deviations: 

1.  At  least  one  hyperplane  giving  the  minimum  sum 
of  absolute  deviations  passes  through  R  of  the  n  points. 

2.  Under  the  assumption  that  no  set  of  R  =  1  observa¬ 
tions  lies  on  one  hyperplane  in  R  dimensions,  if  one 
denotes  n1  as  the  number  of  points  above  a  certain  hyper¬ 
plane  and  n2  as  the  number  of  points  below  this  hyperplane, 
then  the  hyperplane  cannot  minimize  the  sum  of  absolute 
deviations  unless  |n^  -  n2 1  <_  R. 
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3.  If  the  number  of  observations,  n,  is  odd,  then 
any  hyperplane  that  minimizes  the  sum  of  absolute  devia¬ 
tions  passes  through  at  least  one  observation  point. 

4.  For  a  given  set  of  observations,  if  multiple 
hyperplanes  exist  that  minimize  the  sum  of  absolute  devi¬ 
ations,  then  any  convex  combination  of  these  optimal  hyper¬ 
planes  is  also  optimal. 

Charnes,  Cooper,  and  Ferguson  (1955)  demonstrated  the 
applicability  of  simplex  as  a  means  of  minimizing  the  sum 
of  absolute  deviations.  In  their  procedure  they  developed 
a  linear  programming  problem  that  consisted  of  n  con¬ 
straints  where  n  was  equal  to  the  number  of  data  points 
and  that  contained  an  objective  function  that  minimized 
the  sum  of  the  non-negative  deviations  for  each  constraint. 
Wagner  (1959)  modified  the  Charnes  et  al.  formulation  by 
reducing  the  number  of  non-negative  deviations  to  two 
variables  to  solve  the  Chebyshev  criterion:  Minimize  the 
maximum  deviation.  Various  linear  programming  formula¬ 
tions  of  each  minimization  criterion  have  been  addressed 
by  many  authors.  Among  these  are  Barrodale  and  Roberts 
(1973),  Barrodale  and  Young  (1966),  Glahe  and  Hunt  (1970), 
Fisher  (1961),  Sielken  and  Hartley  (1973),  and  Narula  and 
Wellington  (1977) .  Forsythe  (1972)  extended  the  minimum 
sum  concept  to  minimize  the  sum  of  the  deviations  taken 
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to  powers  between  zero  and  infinity  by  the  use  of  gradient 
search  techniques. 

The  preceding  statistically  oriented  models  (i.e., 
least  square  and  least  sum  regression)  were  designed  to 
minimize  the  vertical  distance  between  the  observation 
points  and  the  linear  equation  to  be  estimated.  Using  the 
linear  location  theory  notation,  the  assumption  in  these 
models  was  that  the  minimum  distance  between  the  existing 
facility  i  located  at  a^,  ai2  and  tile  linear  facility, 
x2  =  Cq  +  CjX^,  was  always  equal  to  a^  “  CQ  “  Clail* 

An  exception  to  this  concept  exists  in  the  field  of  sta¬ 
tistical  estimation.  With  the  technique  of  principal  com¬ 
ponents,  the  Euclidean  distance  raised  to  the  second  power 
(p  =  2,  q  =  2)  is  minimized  by  finding  a  linear  composite 
of  the  original  variables  such  that  the  sum  of  the  squared 
perpendicular  distances  from  the  observations  to  the 
linear  composite  are  minimized  (Green,  1978) . 

Summary 

In  summary,  many  solution  techniques  are  available 
to  solve  individual  cases  of  the  generalized,  bivariate, 
linear  location  problem. 

For  regression,  which  minimizes  the  vertical  dis¬ 
tances,  the  least  squares  procedure  is  used  to  minimize 
these  distances  raised  to  the  second  power;  linear 
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to  powers  between  zero  and  infinity  by  the  use  of  gradient 
search  techniques. 

The  preceding  statistically  oriented  models  (i.e., 
least  square  and  least  sum  regression)  were  designed  to 
minimize  the  vertical  distance  between  the  observation 
points  and  the  linear  equation  to  be  estimated.  Using  the 
linear  location  theory  notation,  the  assumption  in  these 
models  was  that  the  minimum  distance  between  the  existing 
facility  i  located  at  a^,  and  the  linear  facility, 

x-  »  C  +  C,x, ,  was  always  equal  to  a.,  -  C  -  C,a.,. 

An  exception  to  this  concept  exists  in  the  field  of  sta¬ 
tistical  estimation.  With  the  technique  of  principal  com¬ 
ponents,  the  Euclidean  distance  raised  to  the  second  power 
(p  =  2,  q  =  2)  is  minimized  by  finding  a  linear  composite 
of  the  original  variables  such  that  the  sum  of  the  squared 
perpendicular  distances  from  the  observations  to  the 
linear  conposite  are  minimized  (Green,  1978) . 

Summary 

In  summary,  many  solution  techniques  are  available 
to  solve  individual  cases  of  the  generalized,  bivariate, 
linear  location  problem. 

For  regression,  which  minimizes  the  vertical  dis¬ 
tances,  the  least  squares  procedure  is  used  to  minimize 
these  distances  raised  to  the  second  power;  linear 
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programming  formulations  are  used  to  minimize  these  dis¬ 
tances  raised  to  the  first  power  and  raised  to  infinity; 
and  gradient  search  procedures  are  used  for  these  dis¬ 
tances  raised  to  any  powers  between  zero  and  infinity. 

For  minimizing  the  Euclidean  distance  (p  =  2),  Wesolowsky's 
procedure  minimizes  p  =  2  distances  raised  to  the  first 
power;  Morris  and  Norback's  modification  of  Wesolowsky's 
procedure,  in  addition  to  being  used  to  minimize  the 
absolute  distances,  can  be  used  for  p  =  2  distances  raised 
to  powers  between  0  and  1;  and  the  principal  components 
method  minimizes  the  square  of  p  =  2  distances. 

Organizational  Flan  of  the  Study 

An  analysis  of  the  generalized,  bivariate,  linear 
location  problem  is  presented  in  this  study. 

Two  methods  for  deriving  the  point  (x^,  x^2)  on  the 
linear  facility  that  minimizes  the  p-norm  distance  taken 
to  the  q  power  to  the  existing  facility  are  discussed 
in  Chapter  II.  In  addition,  solution  properties  and  a 
physical  interpretation  of  the  p-norm  distance  within  the 
generalized,  bivariate,  linear  location  problem  cure  pre¬ 
sented  . 

Contained  in  Chapter  III  is  a  presentation  of  the 

formulation  of  the  modified  distance  function,  . (C  ,  C. )  . 

pi  o  l 

A  modification  of  the  hyperbolic  approximation  technique 
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of  Morris  and  Verdini  (1979)  is  used  to  investigate  the 
convexity  properties  of  J,^(C0#  C1) .  Also,  the  Weber 
formulation  of  the  generalized,  bivariate,  linear  location 
problem  based  on  the  modified  hyperbolic  approximation 
function  along  with  its  convexity  properties  are  discussed. 

Analytical  solutions  of  the  Weber  formulation  of  the 
generalized,  bivariate,  linear  location  problem  for  cer¬ 
tain  values  of  p  and  q  are  presented  in  Chapter  IV.  Also, 
in  Chapter  IV  a  heuristic  solution  for  the  generalized, 
bivariate,  linear  location  problem  for  values  of  p  and  q 
not  analytically  solvable  is  presented. 

Included  in  Chapter  V  are  a  summary  and  recommenda¬ 
tions  for  further  research.  Following  Chapter  V  is  the 
bibliography  which  concludes  this  study. 


CHAPTER  II 
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WEBER  REFORMULATION :  POINT  DETERMINATION 

Consideration  must  be  given  to  procedures  for 
restating  the  Weber  formulation  of  the  generalized,  bivari¬ 
ate,  linear  location  problem  as  presented  in  Chapter  I  in 
a  format  that  can  be  used  to  identify  the  unknowns,  Cq  and 
C^,  of  the  linear  location  facility,  x2  =  Cq  +  C^x.^. 

Locating  a  linear  facility  is  a  two-step  process: 

(1)  determine  the  point  (x^,  xi2)  on  the  linear  facility 
"closest"  to  the  existing  facility  i,  and  (2)  determine 
estimates  for  the  linear  parameters  CQ  and  C^  In  this 
chapter  step  1  is  discussed.  Specifically,  the  point 
given  by  the  coordinates  (x^,  xi2^  minimizes  the 

modified  p-norm  distance  from  the  existing  facility  to 
the  linear  facility  is  determined. 

As  presented  in  Chapter  I,  the  modified  distance 
function,  ilq^(Co,  C^)  ,  is  as  follows: 

o  2  ,a  -  x  ,Plq/p 

Minimize  JT*.  (C  ,  C. )  =  [  I  ait  it  J  (2.1) 

P1  °  1  t=l 

subject  to 

°o  +  cixn  *  xi2' 


p  >.  1,  and 
q  >  0 . 


20 


As  stated,  the  modified  distance  function,  £.^.(C  ,  C, )  , 

pi  o  1 

contains  four  unknowns.  The  obvious  unknowns,  C  and  C. , 

o  1 

are  the  parameters  needed  to  place  the  linear  facility. 

C  and  C.  cannot  be  found,  however,  until  x.,  and  x .  _ , 
o  1  il  iz 

the  points  on  the  linear  facility  that  minimize  the 

p-norm  distance  taken  to  the  q  power  from  the  linear 

facility  to  the  existing  facility  i,  are  determined. 

x..  and  x..~  cannot  be  determined,  however,  until  C  and 
il  i2  o 

are  found. 

In  regression,  which  minimizes  the  sum  of  the  squared 
vertical  distances,  the  values  of  x^  and  x^2  can  be 
stated  as  functions  of  CQ,  C^,  a^,  and  a^2  as  follows: 


x  . .  =  a . .  , 
il  il 


x.0  =  C  +  C,a . . . 
i2  o  1  il 


For  the  Euclidean  (p  =  2)  distances  the  values  of  x^  and 
xi2  minimize  the  modified  distance  function  corres¬ 

pond  to  the  perpendicular  distances  from  the  linear  facil¬ 
ity  to  the  existing  facility  i.  In  general,  however,  the 
values  of  x^1  and  x^2  that  minimize  the  p-norm  distance 
taken  to  the  q  power  are  not  on  the  vertical  or  the  per¬ 
pendicular  and  must  be  analytically  determined. 
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The  first  step  in  solving  the  preceding  circular  con¬ 
dition  is  to  determine  a  general  expression  for  x^1  and 
x^2  as  a  function  of  CQ,  C^,  a^,  and  a^2  for  any  p-norm 
distance  taken  to  the  q  power.  At  this  stage  of  the 
analysis  CQ  and  are  considered  given,  and  the  variables 
to  be  optimized  are  the  coordinates  x^  and  x^2-  In  other 
words,  given  a  line,  x2  =  Co  +  C^x^,  (where  CQ  and  are 
known)  and  a  point  (a^,  ai2)  •  the  coordinates  x^  and 
x^2  on  the  line  that  minimize  the  p-norm  distance  taken 
to  the  q  power  must  be  found.  Mathematically  stated,  the 
problem  is  to  solve  the  following  equation: 


Minimize  (x. ,  x_) 
pi  1  i 


1  2  lait  -  x.  Jp]q/P  (2.2) 

t=l 


subject  to 


Co  +  °lXil  *  x12' 


p  ^  1,  and 
q  >  0. 


Lagrancrietn  Formulation 

One  method  of  finding  the  point  (x^,  x^2)  on  the 
linear  facility,  x2  ■  CQ  +  C^x^,  is  to  solve  the  following 
Lagrangian  formulation: 
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L(X,  X)  =  lji(xil,  xi2) 


+  X (xi2  -  C1xil  -  CQ) . 


(2.3) 


Taking  the  partial  derivatives  with  respect  to  x^,  xi2' 
and  X,  the  following  equations  result: 


=  q|au  -  *illp'1<*’t>1'p/q(-l> 


*  sgn(ai;L  -  x^)  -  XC1# 


(2.4) 


^  41  =  q|a42  -  xf,|P-1(^i)1-p/q(-l) 


axi2 


i2  i2 1  '"pi' 

*  sgn(ai2  -  xi2)  +  X, 


(2.5) 


3L,(X,  X)  =  x  _  c  x  _  c 
3X  Xi2  ^lXil  Lo 


(2.6) 


where  sgn  is  the  signum  function  with  the  following 
properties: 

if  D  >  K,  then  sgn(D  -  K)  =  1, 

if  D  =  K,  then  sgn(D  -  K)  =0,  and 

if  D  <  K,  then  sgn(D  -  K)  =  -1. 

Equations  (2.4),  (2.5),  and  (2.6)  must  each  be  set  equal 

to  zero  and  solved  simultaneously. 

Property  2 . 1 

If  a^  «  x^1  and  a^2  =  x^2,  then  the  first-order 
partial  derivatives  of  L(X,  X)  are  undefined  for  p  >.  q. 
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Proof.  Let  a^  =  x^  and  a^2  =  xj_2  *  Then  the  first- 
order  partial  derivatives  of  L(X,  X)  with  respect  to  x^fc 
(there  t  =  1  or  2)  are  as  follows: 

For  p  =  q: 


=  q|0|p-1[(0)p  +  (0 )p J 1-1  sgn (0)  (-1)  . 

3xit 

For  p  >  q: 

3L(X,  X)  _  j  Q  , p-i  _ -sgn(Q) _ 

3xit  "  '  [  (0)P  +  <0)P]p-^P 

=  q  j  0 | P_1  sgn(0) (1/0) (-1)  . 

Since  [0]°  and  1/0  are  undefined,  3L(X,  X)/3xit  is 
undefined.  The  first-order  partial  derivatives  of  L(X,  X) 
are  undefined  for  p  >.  q. 

Property  2.2 

If  a^t  =  xit  for  t  *  1  or  2,  then  the  first-order 
partial  derivative  of  L(X,  X)  with  respect  to  xit  is 
undefined  for  -1  jc  p  £  1. 

Proof.  Let  a..  -  x.^  for  t  =  1,  2.  Then  the  first- 
-  it  it 

order  partial  derivatives  of  L(X,  X)  with  respect  to  x^t 
are: 

For  p  *  1: 

~Llx!'  ^  °  ql°|1~1(&pi>1~1/q  sgn(0)  (”1)  *  |0|°. 


I 

I 


is 


Again,  since  | 0 | ^  and  1/0  are  undefined,  3L(X,  X)/3x^t 
undefined  for  -1  <  p  <_  1. 

Property  2 . 3 

If  a±1  =  xi;L  and  ai2  *  xi2,  then  the  rth-order  par¬ 
tial  derivatives  of  L(X,  X)  are  undefined  for  p  >  q/r  and 
-r  <  p  <  r. 

Proof.  Using  the  same  arguments  as  those  used  in 
Property  2.2,  the  values  for  the  rth-order  partial  deriva¬ 
tives  of  L(X,  X)  for  the  ranges  of  p  stated  previously 
are  equal  to  0°  or  1/0  and,  therefore,  are  undefined. 

Alternate  Formulation 

The  computational  problems  associated  with  the  signum 
function  and  the  undefined  nature  of  the  Lagrangian  func¬ 
tion,  L(X,  X),  require  that  an  alternate  formulation  of 
the  problem  be  used  to  determine  the  point  (x^,  xi2)  on 
the  linear  facility,  x2  =  CQ  +  cixi'  closest  to  the  exist¬ 
ing  facility  i. 

Again,  recognizing  that  locating  a  linear  facility 
is  actually  a  two-step  process  (i.e.,  the  first  to  deter¬ 
mine  the  point  (x^,  xi2)  on  the  linear  facility  closest 
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to  the  existing  facility  i,  and  the  second  to  determine 
estimates  for  the  linear  parameters  CQ  and  C^) ,  the  problem 
is  divided  into  two  sections.  Step  1  (i.e.,  finding  the 

point  (x^,  as  a  ^unct^-on  CQ>  C]/  aii'  and  ai2 

on  any  linear  facility,  x2  =  Cq  +  cixi'  c^osest  to  the 
existing  facility  i  located  at  (a^,  ai2^  ^etenn;*-ned 
here. 

Since  the  purpose  is  to  determine  the  point 

(x.^,  xi2)  on  a  given  linear  facility,  the  values  of  CQ 

and  can  be  assumed  to  be  known  at  this  stage  of  the 

analysis.  In  theory,  however,  the  process  of  determining 

x ^  and  xi2  would  be  repeated  an  infinite  number  of  times 

for  each  possible  combination  of  values  for  CQ  and  . 

Since  CQ  and  are  assumed  known  at  this  stage  of 

the  analysis,  whether  an  existing  facility  i  lies  on  the 

given  linear  facility  is  known.  If  this  is  the  case,  then 

x . .  =  a . and  x._=C  -  C,a.,  =  a.~  and  the  modified 

il  il  i2  o  1  il  i2 

distance  function  2.q.  equals  zero.  For  an  existing  facil- 

pi 

ity  i  not  on  a  given  linear  facility,  the  Lagrangian 
formulation  for  the  modified  distance  function  as  a  func¬ 
tion  of  xii  and  x^2  is  similar  to  Equation  (2.3) : 

L(X,  X)  =  xi2^  +  X(xi2  ~  Clxil  “  Co)  *  (2.7) 

For  points  on  the  linear  facility,  L(X,  X)  *  0. 
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Substituting  Equation  (2.2)  for  fcq^(x^,  x^2>  •  the 
Lagrangian  formulation  becomes: 

L(X,  X)  =  [ | aiX  -  x±1|P  +  |ai2  -  xi2|p]q/p 

+  X (xi2  -  C1xil  -  CQ) .  (2.8) 


The  computational  problems  associated  with  the 
absolute  value  function  within  the  modified  distance  func¬ 
tion  are  easily  eliminated  by  squaring  and  taking  the 
square  root  of  the  expressions  within  the  absolute  value 
functions.  After  performing  this  operation  and  setting 
the  resulting  relationships  into  a  Lagrangian  formulation, 
the  following  formula  results: 


Minimize  L(X,  X)  =  {[(a.^ 


+ 

+ 


[  (ai2 


X  <xi2 


-  X 

-  X 


\  2 1 p  /  2 
ir  J 


i2 


)  jP/2}q/p 


Clxil  -  Co) 


(2.9) 


Taking  the  partial  derivatives  with  respect  to  each 
variable,  the  following  equations  result: 


+  l(ai2  -  xi2)2]p/2}(q/p)_1 
*  [(ail  '  =‘u>2!<P/2)'1 


*  (aii  ‘  xu>  -  xcr 


(2.10) 
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X,  X)  _ 

3X  "  xi2 


-  C.  x  .  . 
1  il 


(2.12) 


Most  of  the  problems  associated  with  the  original 
Lagrangian  formulation  have  been  eliminated.  The  first- 
order  partial  derivatives  with  respect  to  x^  and  xi2  now 
can  be  computed  since  the  signum  function  has  been 
eliminated,  and  the  first-order  partial  derivatives  are 
defined  for  all  values  of  p  and  q  as  long  as  a^  /  x^ 
and  a^2  #  xi2*  T^^-s  nonequality  condition  for  both  vari¬ 
ables  has  been  eliminated.  One  problem,  however,  still 

remains.  If  just  one  a..  =  x .  (i.e.,  either  a..  =  x.. 

it  it  il  il 

or  a^2  =  xi2^  '  then  3L(X,  X) /3x^t  is  undefined  for  p  <_  2. 

Ignoring  this  property  for  the  moment  and  recognizing 
that  the  following  derivation  may  not  be  valid  for  p  <_  2, 
the  values  of  x^1  and  x^2  can  be  solved  by  setting  Equa¬ 
tions  (2.10),  (2.11),  and  (2.12)  equal  to  zero  and  solving 
simultaneously.  Performing  these  steps  and  eliminating 
the  following  equation  results: 
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Uu  -  x^X.a^  -  xu,2]<P/2>-1 


=  -C1(ai2  "  CQ  -  C1xil) 


t'ai2  -  Cl* il  -  C0)2l(p/2>-1. 


(2.13) 


Squaring  both  sides  of  the  equation  to  ensure  that  is 
positive  and  taking  the  (p-l)th  root,  the  following  formula 
is  obtained: 


,ail  -  aii>2  ■ 


(ai2  “  Co  Clxil)  ‘ 


(2.14) 


Expanding  and  collecting  terms  and  then  using  the  quad¬ 
ratic  formula,  the  following  values  of  x^  are  obtained: 

,ail  HC12,1/2]1/(p-1>  iai2  -  cl 

1  -  tci2JP/ ‘p-1’ 


{1  -  [Cx] ((C12)1/2]1/(P"1)} , 


(2.15a) 


an  -  [,Cl2>l/2,l/.P-»t.<2 

1  -  [c^jP/'p-11 

{1  +  [C11[(C12)1/2]1/(P_1)}. 


Cl 

_o_  } 


(2.15b) 


Substituting  the  values  of  x^  into  Equation  (2.12),  the 
following  values  for  x^2  are  obtained: 
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C  +  C  a  —  FC 

co  +  ciaii  LCi  3 _ til 

1  -  [C1V/,P-:U 

ICintC12)1/2ll/,p'1’ [ai2  -  Cc 
l  - 


Clail' 


(2.16a) 


C  +  C  a  —  rp  ^iP/  {P~^-)a 
co  +  ciaii  LCi  1 _ til 

i  -  [c12)p/(p_1) 

fC1]  [(Cj2)172]37^1*  [ai2  -  CQ  - 

i  -  Ic12]p/<p-1> 


Clail] 


(2.16b) 


The  pair  of  solutions  for  and  results  from 
the  squaring  of  in  the  solution  procedure.  This  tech¬ 
nique  is  necessary  since  the  slope  (C^)  may  be  negative  and 
the  (p-l)th  root  may  be  noninteger.  Resolving  the  con¬ 
flict  of  two  solutions  for  x^  and  xi2  is  accomplished 
by  substituting  each  pair  of  values  into  Equation  (2.13). 
Performing  this  action  for  the  first  values  of  x^  and 
xi2#  the  following  equation  results: 


2)1/2J  =  -cr 


(2.17) 


The  term  on  the  left  side  of  the  equal  sign  is  always  less 
than  or  equal  to  zero,  and  the  only  way  for  the  equation 
to  be  equal  is  for  C1  _>  0.  Thus,  the  first  values  of  x^ 
and  x^2  are  valid  solutions  if  and  only  if  >  0. 


1 
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Likewise,  performing  the  same  action  with  the  second 
values  of  x^  and  x^2,  the  following  equation  results: 


t (CL2) 1/2]  =  -C1.  (2.18) 

In  this  case  the  term  on  the  left  side  of  the  equation 
sign  is  always  greater  than  or  equal  to  zero.  Thus, 
must  be  less  than  or  equal  to  zero,  and  the  second  values 

of  and  x^2  are  valid  solutions  if  and  only  if  c^  <  0. 

/ 

Solution  Properties  of  x^  and  x^2 

The  following  properties  and  associated  proofs  sub¬ 
stantiate  the  assertion  that  the  point  (x^,  x^2)  mini¬ 
mizes  the  modified  p-norm  distance  from  the  existing 
facility  i  to  the  linear  facility. 


(2.19 
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+  (p  -  1){  Ua±1  -  xu)2]p/2 

*  [<ai2  '  Co  -  Cl*il,2lP/2} 

*  {[(an  -  xil)2]  ^/2)~1 

+  Cl2[(ai2  '  Co  "  c1*il)23(P/2,"1}l 


(2.21) 


Notice  that  the  portion  of  Equation  (2.21)  that  is  under¬ 
lined  is  equal  to  zero  for  solutions  derived,  and  all  of 
the  remaining  terms  are  strictly  positive  except  those 
containing  q  and  (p-1) .  Thus,  the  second  derivative  of 
L(x^)  is  strictly  positive  as  long  as  q  and  (p-1)  are 
greater  than  or  equal  to  zero.  The  second  derivative  is 
also  strictly  positive  if  q  >  0  and  p  >_  1. 

Property  2.5 

The  point  (x^,  x^2)  on  linear  facility, 
x2  =  CQ  +  C^x^,  that  minimizes  the  modified  distance 
function  .  (X)  is  independent  of  q  as  observed  from 
Equations  (2.15a),  (2.15b),  (2.16a),  and  (2.16b). 

Property  2.6 

If  the  existing  facility  i  is  on  the  linear  facility, 

x2  "  Co  +  C1X1?  then  Xil  =  ail'  and  Xi2  =  ai2* 

Proof.  (Since  the  proofs  for  x^  *  a^1  for  Equa¬ 
tions  (2.15a)  and  (2.15b)  and  the  proofs  for  x^2  *  ai2 


for  Equations  (2.16a)  and  (2.16b)  are  similar,  only 


=  a^1  for  Equation  (2.15a)  is  proven  here.)  If  the 
existing  facility  i  is  on  the  linear  facility,  then 
a^j  =  CQ  +  C^a^.  Substituting  this  relationship  into 
Equation  (2.15a),  the  following  equation  results: 


an  +  [<c12)1/2]1/(p-11  ;c0  +  ciail 

2.p/(p-l) 


kil 


Cl 

_o_} 


i  -  ic^r 


*  (1  -  [  (C12) 1/2]  -17  (P“1>  >  . 


(2.22) 


Since  must  be  greater  than  zero  for  Equation  (2.15a) 
to  be  a  solution  to  Equation  (2.13),  the  following  opera¬ 
tion  can  be  performed: 


Xil 


au{l  +  [(C12)1/2]p/(p"I)Hl  -  [(C12)1/2]p/(p~1)} 


1  -  (c1"2ip7Tprrr 


*  ail* 


The  following  properties  eliminate  the  problem  associ¬ 
ated  with  the  first-order  partial  derivatives  of  the  modi¬ 
fied  distance  function. 


Property  2.7 

The  slope  of  the  line  connecting  the  existing  facility 
i  to  the  point  (x^,  x^2)  on  the  linear  facility, 
x2  *  CQ  +  CjX^,  that  minimizes  the  p-norm  distance  taken 


* 
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to  the  q  power  is  equal  to  -1/ [  (C^2) 1//2]  ly^p  ^  if  C1  >  0 
and  equal  to  1/ l (Cj2) 1/2] (p-1)  if  CL  <  0. 

Proof.  As  previously  stated,  the  coordinates  of  the 
point  on  the  linear  facility  that  minimize  the  p-norm 
distance  taken  to  the  q  power  from  the  linear  facility, 
x2  =  CQ  +  C-^x^,  to  the  existing  facility  i  when  >_  0  are 
given  by  Equations  (2.15a)  and  (2.16a).  Substituting 
these  values  into  the  definition  of  the  slope,  the  follow¬ 
ing  equation  results: 


Slope  = 


AX, 

d 

AX. 


ai2  ~  xi2 _ zl _ 

ail  “  Xil  t(C12)1/2]1/p_1’ 


For  C-l  £  0  the  optimal  values  for  and  xi2  are 
given  by  Equation  (2.15a)  and  Equation  (2.16b),  respec¬ 
tively.  Substituting  these  values  into  the  definition  of 
the  slope,  the  following  equation  results: 


ax2 

slope  "  siq; 


ai2  ~  xi2  _ 1 _ 

aii  "  xii  '  [(c12i1/2]1/P-1‘ 


If  the  existing  facility  i  is  not  on  the  linear 
facility,  x2  *  CQ  +  CjX^  and  /  0;  then  x^  =  a^  if 
and  only  if  p  equals  the  limit  as  p  approaches  1  from  the 
negative  direction. 

If  the  existing  facility  i  is  not  on  the  linear 
facility,  x2  85  CQ  +  cix]/  then  the  only  way  for  x^  =  a^1 


if  slope  of  the  line  connecting  the  two  points  is 
vertical  (i.e.,  AXj/AX^^  approaches  infinity).  Since  the 
only  way  the  slope  can  be  equal  to  infinity  is  when 
(1-p)  *  0  (i.e./  the  limit  of  (1-p)  as  p  approaches  1 

from  the  negative  direction) ,  then  the  limit  of  the  slope 
as  p  approaches  1  from  the  negative  direction  is  as  fol¬ 
lows  : 

lim_  ( - — - j— ) 

1/P 

=  ^‘-(-[(C12)1/2]1/P  ‘  P/P  = 

Property  2.8 

If  the  existing  facility  i  is  not  on  the  linear 
facility,  x2  =  CQ  +  Cj^,  and  Cj  *  0;  then  x±2  =  a±2  if 
and  only  if  p  equals  the  limit  as  p  approaches  1  from  the 
positive  direction. 

Proof .  If  the  existing  facility  i  is  not  on  the 
linear  facility  and  ^  ii  0,  then  the  only  way  xi2  =  ai2 
is  when  the  slope  of  the  line  connecting  the  two  points 
is  horizontal  (i.e.,  AXj/AXj^  *  0).  Taking  the  limit  as 
P  approaches  1  from  the  positive  direction,  the  following 


formula  results: 
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lim  , 
p-l+' 


nc12)1/2ip"1 


) 


[ (c12)1/2]p/p  ‘  1/p 


Thus,  the  property  is  proven. 

As  can  now  be  seen  using  Property  2.7,  Equation  (2.11) 
is  defined  for  all  values  of  p  and  q  with  the  exception  of 
p  =  1.  Likewise,  as  can  be  seen  using  Property  2.8, 
Equation  (2.12)  is  defined  for  all  values  of  p  and  q 
except  p  =  1.  The  first  derivatives  are  undefined  at 
p  =  1  because  the  Lagrangian  function  is  discontinuous  at 
p  =  1.  The  point  (x^,  x^)  that  minimizes  the  modified  p 
distance  from  the  existing  facilities  to  the  linear  facil¬ 
ity  for  p  =  1  is  as  follows: 

Limit  as  p  approaches  1  from  the  negative  direction: 


xi2  ”  ai2 ’ 

Thus,  the  derivation  of  the  point  (x^,  x^2) 

for  all  values  of  p  and  q,  and  the  point  (x^,  x^) 
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minimizes  the  p-norm  distance  taken  to  the  q  power  with 
q  >_  0  and  p  >  1  or  when  q  <  0  and  p  <  1. 

Interpretation  of  x^,  x^2 

The  results  of  Property  2.7  indicate  a  physical 
interpretation  for  p  within  the  modified  distance  func¬ 
tion  Jlpi*  Working  only  with  the  case  where  >  0  (even 
though  the  conclusions  are  generalizable  to  <  0)  ,  the 
slope  of  the  line  connecting  the  existing  facility  to  the 
linear  facility  is  equal  to  -1/  [  (C^)  ^2]  ^  .  That 

the  slope  of  this  line  is  either  vertical  or  horizontal 
as  p  approaches  1  from  the  negative  or  positive  direction 
has  been  shown.  That  the  slope  of  this  same  line  is  equal 
to  -1  as  p  approaches  infinity  can  also  be  shown.  For 
any  other  value  of  p  the  slope  of  the  line  minimizing  the 
modified  p-norm  distance  is  a  function  of  the  relation¬ 
ship  between  the  values  of  p  and  C^.  If  p  =  2,  for 

2  1/2 

example,  the  slope  is  equal  to  -1/[(C^  )  ]  which  is  the 

perpendicular  distance  when  >  0.  As  p  gets  larger, 
the  absolute  value  of  the  slope  of  the  line  connecting 
(ail»  »i2)  to  the  linear  facility  for  a  given  value  of 
C1  gets  larger. 

Since  the  p-norm  distance  taken  to  the  q  power  is 
minimized  in  the  generalized,  bivariate,  linear  location 
problem,  the  value  of  p  determines  which  distance  will  be 
minimized.  When  p  approaches  1  (i.e.,  the  rectalinear 


%mm  ■■  t 
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case),  the  minimized  distance  is  in  the  direction  of  only 
one  of  the  variables.  To  use  regression  terminology, 
minimizing  the  error  in  only  one  direction  assumes  that 
one  variable  is  fixed  and  that  the  other  variable  is 
random  and  accountable  for  the  total  error  (i.e.,  indepen¬ 
dent  and  dependent  variables) .  When  p  is  not  equal  to  1, 
the  error  is  divided  between  the  two  variables  according 
to  -1/  [  (C.^)  and  the  terms  "independent"  and 

"dependent"  variables  are  no  longer  appropriate.  Thus,  p 
is  the  parameter  within  the  generalized,  bivariate,  linear 
location  function  that  determines  the  relative  weight 
between  the  two  variables  with  respect  to  the  error  term. 


CHAPTER  III 


WEBER  REFORMULATION:  MODIFIED 
DISTANCE  FUNCTION 


As  previously  stated,  locating  a  linear  facility  is 
a  two-step  process  involving  determination  of  (1)  the 
point  (x^,  xi2)  on  the  ^^-near  facility  "closest"  to  the 
existing  facility  i,  and  (2)  estimates  for  the  linear 
parameters  CQ  and  C. .  In  Chapter  II  a  method  of  deter¬ 
mining  the  point  (x^,  xi2)  on  a  line,  x2  =  Cq  +  C.^, 
that  minimizes  the  modified  p-norzn  distance  to  the  exist¬ 


ing  facility  i  located  at  the  point  (a^, 
trated . 


ai2)  was  illus 


Before  estimates  for  the  linear  parameters  CQ  and  C^ 
can  be  determined,  more  detailed  consideration  must  be 


given  to  the  Weber  formulation  of  the  generalized,  bivari¬ 
ate,  linear  location  problem.  In  this  chapter,  therefore, 

the  modified  distance  function  (C  ,  C, )  associated 

pi  o  l 

with  the  p-norm  distance  taken  to  the  q  power  from  the 
existing  facility  i  to  the  point  (x^,  x^2),  on  the  ^inear 
facility  is  developed.  The  convexity  or  lack  of  convexity 
of  is  also  investigated.  Finally,  the  Weber 

formulation  of  the  linear  location  problem  is  presented. 
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Formulation  of  (CQ,  C^) 

The  values  of  x^  and  x^  developed  in  Chapter  II  can 
now  be  used  to  formulate  the  modified  distance  formula 
^.<co,  C^)  .  Inserting  these  values  into  Equation  (2.1), 
the  following  two  formulations  corresponding  to  the  pairs 
of  values  for  x^  and  x^  result: 


*pi(Co'  Cl»  • 


[U  -  [ClH(C12)1/2l1/(p'1,}2)q/2 

tu  -  (c12ip/(p‘1,}2)q/2 
(1  +  [(c12)1/2]p/<P"l)}'5/p 

[(ai2  -  co  -  Vii,2iq/2' 


(3.1a) 


[{1+  [C11 [(ci2)1/2]1/(P"1,}2]P/2 
[{l  -  [c12]p/(p_1)}2]p/2 

*  a  +  [(c12)1/2]p/(p-1)}q/p 


*  “ai2  -  Co  -  C1ail,2>q/2' 


(3.1b) 


The  solution  in  Equation  (3.1a)  resulted  from  using 
the  first  solutions  for  x^  and  x^  (i.e.»  Equations 
(2.15a)  and  (2.16a)),  and  the  second  solution  in  Equa¬ 
tion  (3.1b)  .^suited  from  use  of  Equations  (2.15b)  and 
(2.16b).  The  only  difference  between  Equations  (3.1a) 
and  (3.1b)  is  found  in  the  following  terms: 


[{1  -  [c1inc12>1/2]1/<p-1>}2iq/2>, 

[U  +  [Cj  [(C12)1/2]1/(p_1)}2]q/2. 


(3.2a) 


(3.2b) 


Remembering  that  the  first  pair  of  solutions  is  valid  only 
when  >_  0  and  that  the  second  pair  of  solutions  is  valid 
only  when  <_  0,  the  following  operations  can  be  per¬ 
formed: 

[{1-  IcJ  [(C12)1/231/(p'1)}2]q/2,  (3.3a) 

[{1-  | C1 i [(C12)1/2]1/(p"1)}2]q/2.  (3.3b) 

Thus,  the  two  terms  are  identical,  and  the  following 

single  formulation  for  (C  ,  C. )  results: 

pi  o  1 

a  t<ai2  -  C  -  Cxa  )2]q/2 

pi*  o'  V  (1  +  [(c  2)l/2]p/(p-l)]q(p-l)/p* 

A  simplified  definitional  formula  for  ilq.  is  as  follows: 

pi 


*pi<Co'  Cl>  " 


Iai2  ~  Co  ~  Claill 

i  +  |C  |P/(p-»|q(p-i)/p 


(3.5) 


Convexity  Properties  of  Lq^ 

To  investigate  the  convexity  of  the  modified  distance 

function  4p^(CQ,  C^)  the  first-  and  second-order  partial 

derivatives  of  2,q .  (C  ,  C. )  with  respect  to  C  and  C,  must 
pi  o  l  o  l 
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be  taken  (Hillier  &  Lieberman,  1980).  Using  either  the 
calculation  Formula  (3.4)  or  the  Definitional  Formula 
(3.5),  similar  difficulties  to  those  encountered  in 
Equation  (2.5)  are  encountered  with  respect  to  undefined 
rth  partial  derivatives  if  q  <  r  and  a^2  =  CQ  + 

To  remove  these  problems,  a  smoothing  constant  m  similar 
to  the  concept  used  by  Morris  and  Verdini  (1979)  for  the 
point  location  problem  can  be  employed.  Using  the  hyper¬ 
bolic  smoothing  constant  m,  the  following  equation  is  an 
approximation  for  &^(CQ,  C^)  : 


lpi«V  ci>  - 


[Vj2  +  m]q/2 
Gq(P-i) /p 


(3.6) 


One  definition  and  several  theorems  from  Protter  and 
Morrey  (1964)  are  helpful  with  respect  to  the  proofs  of 
properties  within  this  section. 


Definition  3.1 

A  function  f  is  said  to  be  continuous  at  (b,  d)  if 
and  only  if : 

(a)  f  is  defined  at  (b,  d) ,  and 

(b)  lim  f(x,  y)  =  f(b,  d) . 

(x,  y)-Mb,  d) 

Theorem  3.1 

If  b  is  a  constant  and  f(x)  =  b  for  all  values  of  x 
then  for  any  number  d 

lim  f (x)  =  b. 
x-*d 


Theorem  3.2 

If  m  and  b  are  constants  and  f(x)  =  mx  +  b  for  all 
values  of  x,  then  for  any  number  d 

lim  f(x)  =  md  +  b. 
x-*d 

Theorem  3.3 

If  f  and  q  are  two  functions  with 

lim  f(x)  =  M.  ,  lim  q(X)  =  M_, 
x-*d  1  x-^d  * 


lira  f (x) 
x-+d  g(x) 


M1 

jj-  where  M2 


*  0 


then 
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and 


lim  [f(x)  *  g(x)]  =  M,  *  M. 
x+d  1 


Theorem  3.4 


If  f  and  g  are  two  functions  with 

lim  f(x)  =  M.  lim  g(x)  =  M_, 

xr*d  1  x+d  * 


then 


lim  tf(x)  +  g(x)3  =  M,  +  M, 
x-*-d  A 


Theorem  3 . 5 

If  n  is  any  positive  integer,  and  if  f  is  such  that 

lim  f(x)  =  M; 
x+d 

then 

lim  [f(x)]1/n  =  M1/n. 
x-*d 


Property  3.1 


The  limit  of  Lq^(CQ,  C^)  as  m  approaches  zero  is 


*pi<Co'  Cl>- 


Proof .  Using  theorem  3.1,  the  following  can  be 


lim  V.2  *  v  2  and  lim  Gq(p"1)/p  -  Gq(p_1,/p  -  L*. 
nr*0  1  1  m+0  1 


stated : 
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Using  theorem  3.2,  the  following  can  be  stated: 


lim  [V.2  +  m]  =  V.2  =  L*. 
nr*0  1  1  z 


Finally,  using  theorems  3.3  and  3.5  and  the  knowledge  that 
q/2  may  be  represented  by  a  rational  number  without 
severely  limiting  the  domain  of  q,  then 


u*  Cl> 


m+0 


lim  KV*2 
irt+0  1  L* 


> 


tL2] 

L1 


q/2 


=  £pi(co,  C^).  (Q.E.D.) 


The  first-  and  second-order  partial  derivatives  of 

Lq.(C  ,  C, )  with  respect  to  C  and  C.  (i.e.,  Equation 
pi  o  1  o  1 

(3.6))  are  as  follows: 


3Lpi(co-  V  .  .  -q  ,  [V  2  +  ml  (<2/2)-1v  ,  (3-7» 

3C  Lq(p-l)/pJ  1  i  J  i 


^oi^o'  Cl)  _  ,  a  2  , (q/2) -2 

172 - ofTi P“ i")  7p  fVi  +  ral 


*  [V^tq-l)  +  m]. 


(3.8) 


.Ww'< 
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3Lgi(Co'  Cl> 


[q (p-1) /p] +11  lvi 


]  [V,2  +  in]  (<I/2)~1 


{GVii  +  cil(ci2>1/2i(2'P)/(P'1) 


*  (V^  +  m)  1)  , 


(3.9) 


9  Si(Co'  Cl>  _  ,q[Vi 
n  —  t  nr 


2  +  m]«*/2>-2 


T-1-HG2ail2[Vi2(q-l)  +  m] 


[q (p-1) /p] 


+  2qCailViG[ (C12)1/21 (2“P)/(P~1) [Vi2  +  m] 

[<C  2)l/2](2-p)/(p-l)[v  2  +  m]2 

+  — 1 - <P^ii - 

*  [ (q  +  1) (p  ~  1) r(C12)1/2]p/(p~1)  -  1), 


(3.10) 


^pi^o'  Cl)  q[Vi2  +  »l(q/2)"2 
_ P1  ° _ ±_  *  r  _ \ 

<■ r  t  tr^—t  \  /v.1  J.1  1 


3C  3C. 
o  1 


*  (  ^(p-D/pl+l  )(S«lltVl  O-11  +  ml 
+  qC1[(Cl2)1/2]2'p/P‘1V1(Vi2  +  ml).  (3.11) 


Property  3.2 


Lq^ (Cq/  C^)  is  continuous  for  q  >.  0  and  p  >  1, 

Proof.  Let  f^(CQ,  C^)  *  Lpi^Cof  C]^  '  and  consider 

f . ( C,  C, )  at  any  point  c  e  C  and  c.  e  C,  where  c  and 
1  O  X  OO  1  1  o 

are  real  numbers.  By  inspection,  requirement  (a)  of 


.  '•  -  »*i  '  ri,  ' 


Definition  3.1  is  satisfied  since  f^(cQ,  c^)  =  (V2  +  m]<^2/ 
^(P'D/p,  i  =  1,  .  ..,  n  are  defined  for  all  cq  s  Cq  and 
c.^  e  C1  and  q  >_  0  and  p  >  1.  Requirement  (b)  is  satisfied 
because 


(b,  d)-McQ,  c^)  i 


f ± (b ,  d)  =  fi(cQ,  c ^ 


for  all  i  *  1,  . ..,  n  and  c  e  C  and  c.  e  C.  .  The  proof 

o  o  11 

for  the  preceding  equation  is  similar  to  the  proof  in 
Property  3.1.  Therefore,  the  function  is  continuous  for 
q  ^  0  and  p  >  1. 

The  case  of  p  =  1  is  not  covered  by  the  preceding  and 
subsequent  properties  because  f^(Co,  ^  is  undefined  at 
p  *  1.  The  situation  where  p  =  1  is  addressed  separately 
with  properties  of  the  Weber  formulation. 


Property  3.3 

The  first-  and  second-order  partial  derivatives  of 

L^. (C  ,  C. )  with  respect  to  C  and  C,  are  continuous  for 
pi  o  l  o  l 

q  >.  0  and  p  >  1. 

Proof.  Let  fi0<C0,  C1)-SLji(C0,  C^/3^, 

-  SLpi<Co'  Cl»/8C1'  f"io  *  ^Lpt(co'  cl)/3Co2'  fil(Co'  cl> 

-  32L’i,(C0.  C^/SCj2,  and  fjl0(C0,  C1>-32^1(Co,C1I/3C1C0 

t  •  n 

and  consider  fiQ(C0,  C^,  fil(Cc>,  ci>,  fio(Co'  ci)  ' 

It  If 

fil*Co'  Cl*  and  filo*Co'  Cl*  at  any  Point  c0  e  Co  and 
C1  e  C2  where  c0  3^1^  c1  are  real  numbers.  Again,  by 


inspection  fio<cQ,  aj  ,  cj  ,  fio(cQ,  , 

n  n 

f^1(cQ,  c^),  and  f .1n(cn,  c.,)  satisfy  requirement  (a)  of 


ilo'  o'  ~l 


Definition  3.1  since  they  are  defined  for  all  c  e  C  and 

o  o 

gx  e  C1#  i  =  1,  ...»  n  and  q  >  0  and  p  >  1.  Requirement 
(b)  is  satisfied  because 


(b,  d)-*-  (cQ,  cx) 


d)  = 


3Lpi*co'  cl) 


(b,  d)-*-(cQ,  c^) 


fn(b,  d)  = 


_  3LDi(co-  el> 


(b,  d)-McQ,  Cj^) 


fio<b'  d> 


9  Lgi«V  cl> 


(b,  d)-*- (cq,  c^ 


f±1(b,  d)  = 


°i’ 


«  32L^.(co,  c,) 

lim  fiio(b'  d> - 1 tic - ~ 

(b,  d)+(c  ,  c±)  110  9Cl9Co 


for  all  i  *  1,  n  and  c  e  C  and  c,  e  C, .  The  proofs 

o  o  11 

for  the  preceding  equations  are  similar  to  the  proof  in 
Property  3.1.  Therefore,  the  property  is  proven. 

Property  3.4 

Lpi (Cq r  C^)  is  strictly  convex  with  respect  to  C0  for 
q  >.1  and  m  >  0  and  strictly  concave  for  0  <  q  <  1  and  the 
limit  as  m  approaches  but  does  not  equal  zero. 


•  -*«  **  Mferto-.v- 


1  V  ,-s '  /Jn*J X 
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Proof .  Since  L^(Cq,  C^)  has  continuous  second 

derivatives  for  all  points  c  e  C  where  c  is  a  real 

o  o  o 

number,  it  is  sufficient  to  show  that  the  second-order 

partial  derivative  with  respect  to  CQ  is  strictly  positive 

for  q  >  1  and  negative  for  0  <  q  <  1  where  the  limit  of 

m  approaches  but  does  not  equal  zero.  From  equation 

(3.8)  the  second-order  partial  derivative  of  Lg.(C  ,  C, ) 

pi  o  1 

with  respect  to  CQ  is  given  by 


32^i(Co-  Cl> 

3C  2 


,q(p-l)/plvi 


[V.2  +  m]  (q/2) ”2 


*  [V±4 (q-1)  +  m] . 


Since  G  and  m  are  always  positive,  3  Lq. (C  ,  C,)/3C  is 

pi  o  1  o 

strictly  positive  for  q  ^  1.  If  0  <  q  <  1,  32Lp^  (CQ,  C^)/ 
2 

3Cq  is  strictly  negative  as  m  approaches  but  does  not 
equal  zero . 


Properties  of  the  Weber  Formulation 
Substituting  Equation  (3.5)  into  the  Weber  formula¬ 
tion  produces  the  following  generalized,  bivariate,  linear 
location  problem: 


1 :"ilai2  -  =0  -  Vul* 

CVCo'  Cl>  ’  (1  +  |Ci|P/(p-l)j‘!(p-l)/p' 


(3.12) 


Again  using  the  hyperbolic  approximation  function  to 
ensure  defined  derivatives  and  to  eliminate  computational 
difficulties  associated  with  the  absolute  value  functions. 


the  following  approximation  of  CJZ.p(CQ,  C1)  results  from 
substituting  Equation  (3.6)  into  the  Weber  formulation: 


clp(co'  cl>  = 


[1  ♦  [,c2,1/2,P/(P-l)1qlP-:1)/p- 


(3.13) 


Earlier  the  statement  was  made  that  the  previous 
properties  in  this  chapter  were  applicable  only  when  p  >  1. 
The  following  properties  relieve  this  restriction  for 

p  -  1. 


Property  3.5 

2 

The  limit  of  CJlp(Co,  C^)  as  p  approaches  1  from  the 
negative  direction  is  the  regression  of  x2  on  x^  taken  to 
the  q  power. 

Proof .  The  limit  of  J as  p  approaches  1 

from  the  negative  direction  is  as  follows: 


lim  (1  ♦  |<v  |p/(p_1)]  -  [1  +  0] 0  -  1.  (3.14) 

P-1-  1 

Thus,  the  Weber  formulation,  when  p  approaches  1  from  the 
negative  direction,  is  as  follows: 
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C£’(C  , 
P  o 


Cl>  = 


n 

z  w. 

i=l  1 


ai2 


-  C  - 

o 


Vill 


(3.15) 


When  Wi  *  1  for  i  =  1,  . ..,  n  and  q  =  2,  the  generalized, 
bivariate,  linear  location  problem  reduces  to  the  least 
squares  regression  of  x2  on  x^. 

Property  3.6 

The  limit  of  C£^(Cq,  C^)  as  p  approaches  1  from  the 
positive  direction  is  the  regression  of  and  x2  taken 
to  the  q  power. 

Proof .  From  Equation  (3.15)  C£^(Cq,  )  is  equal  to 

the  following  equation: 


C£^(C  ,  )  = 

p  o  1 


n 

Z  W. 
i=l  1 


ai2 


”  Co  " 


Clail> 


[l  +  |c1|p/(p~1)]q(p"1)/p 


Dividing  the  numerator  and  denominator  by  l/|C.jJq, 

Ca5(C„,  C. )  becomes 
p  o  1 


c*plco-  Cj) 

Taking  the  limit  of 
positive  direction. 


n 

Z  W. | (a 
i»l  1 


i2 


-  co  -  clail,/c 


Cip(Co#  Cx)  as  p  approaches  1  from  the 
the  following  result  is  obtained: 
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lim  [CS,^(C 
p-l+  p 


o' 


^1Wil<ai2-Co-Clail)/Cllq- 


(3.16a) 


Rearranging  the  linear  constraint,  x2  =  CQ  +  cix]/  to 

I  ( 

x^  =  (x2  -  CQ) /C1  and  defining  =  1/C1  and  CQ  =  ~CQ/Ci' 
the  limit  of  C£^(CQ,  C.^)  as  p  approaches  one  from  the 
positive  direction  may  be  stated  as 


Equation  (3.16b)  is  just  the  statement  of  the  objective 
function  for  regression  x^^  on  x2*  As  can  be  seen,  regres¬ 
sion  of  x^  on  x2  may  be  solved  indirectly  as  a  special 
case  of  the  generalized,  bivariate,  linear  location 
problem.  The  preceding  result  is  supported  by  Property 
2.7  where  the  statement  was  made  that  the  slope  of  the 
line  connecting  the  existing  facility  i  to  the  point  on 
the  linear  facility  that  minimizes  the  p-norm  distance 
taken  to  the  q  power  is  horizontal  (i.e.,  equal  to  zero) 
when  the  limit  is  taken  as  p  approaches  1  from  the 
positive  direction. 
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Property  3.7 

CLp  ( CQ t  is  strictly  convex  with  respect  to  Cq  for 

q  ^  1  and  m  >  0  and  strictly  concave  for  0  <  q  <  1  and  the 
limit  as  m  approaches  but  does  not  equal  zero. 

Proof.  This  property  follows  since  CLp(CQ,  is 

the  sum  of  strictly  convex  terms  when  q  _>  1  and  m  >  0  and 
strictly  concave  terms  when  0  <  q  <  1  and  the  limit  as  m 
approaches  but  does  not  equal  zero. 

No  general  statement  can  be  made  at  this  point  of  the 
analysis  concerning  the  overall  convexity  or  lack  of  con¬ 
vexity  of  CL^(Cq/  C^)  since  this  property  changes  with 
different  values  of  p  and  q.  By  Property  3.7  CLq(C  ,  C,) 
is  not  convex  with  respect  to  CQ  for  values  of  0  <  q  <  1. 
This  is  a  sufficient  condition  for  stating  that  CL^(Cq,  C^) 
is  not  convex  within  the  range  0  <  q  <  1  (Hillier  & 
Lieberman,  1980).  Yet,  by  Properties  3.5  and  3.6 
CLp(CQ,  C^)  reduces  to  simple  regression  when  q  =  2  and 
the  limit  is  taken  as  p  approaches  1  in  either  the  plus 
or  minus  direction  and  simple  regression  is  strictly  con¬ 
vex  (Neter  &  Wasserman,  1974).  This  inconsistency 
requires  that  convexity  be  investigated  for  each  indi¬ 
vidual  combination  of  p  and  q. 


Summary 

Detailed  in  this  chapter  is  the  Weber  formulation  of 
the  generalised,  bivariate,  linear  location  problem.  The 
problems  associated  with  undefined  derivatives  require  the 
addition  of  a  hyperbolic  smoothing  constant.  The  modified 
Weber  formulation  of  the  bivariate ,  linear  location  problem 


is  not  necessarily  concave  or  convex. 


CHAPTER  r V 


SOLUTIONS  TO  THE  GENERALIZED,  BIVARIATE, 

LINEAR  LOCATION  PROBLEM 

The  complex  structure  of  the  Definitional  Formula, 
Equation  (3.12),  and  the  Weber  Reformulation  Equation 
(3.13),  along  with  the  associated  nonconvexity  prohibits 
a  universal  analytical  technique  for  solving  the  general¬ 
ized,  bivariate,  linear  location  problem.  Yet,  the  situa¬ 
tion  is  not  hopeless  for  those  requiring  analytical 
solutions.  Just  as  illustrated  in  Chapter  III,  linear 
regression  (i.e.,  q  *  2  and  p  approaches  1  from  either  the 
positive  or  the  negative  direction)  is  one  special  case 
of  the  generalized,  bivariate,  linear  location  problem. 
Other  special  combinations  of  p's  and  q's  likewise,  reduce 
to  common  forms  of  analysis  that  are  currently  employed. 

In  this  chapter  the  relationship  between  the  general¬ 
ized,  bivariate,  linear  location  problem  and  current  forms 
of  analysis  is  illustrated.  In  addition,  certain  unique 
combinations  of  p's  and  q's  that  represent  models  not 
presently  employed  are  solved  by  various  analytical  or 
numerical  techniques.  Finally,  a  heuristic  technique  is 


presented  for  a  general  solution  of  the  generalized, 
bivariate,  linear  location  problem. 


The  Limit  as  p  Approaches  1  Family 
If,  in  a  manner  similar  to  that  used  in  the  proof  of 
Property  3.5,  the  limit  is  taken  as  p  approaches  1  from 
the  negative  direction,  the  definitional  formulation  of 
the  generalized,  bivariate,  linear  location  problem  Equa¬ 
tion  (3.12)  becomes  the  following  formula: 


C£?-(C 


Cl)  = 


n 

E  W.|a 
i=l  1  1 


i2 


'  Co  ' 


Claill 


(4.1) 


The  distance  being  considered  in  this  family  of  models  is 
the  vertical  distance  from  the  existing  facilities  to  the 
linear  facility  to  be  located  taken  to  the  power  q. 

As  shown  in  Property  3.5  when  q  =  2,  Equation  (4.1) 
becomes 


Ctl— (co' 


n 

c. )  =  E  w.  (a 
1  i=l  1 


i2 


-  Co  - 


Clail> 


(4.2) 


which  is  simply  the  regression  of  x2  on  x^  The  solution 

for  this  well-known  model  is  found  by  using  the  classical 

unconstrained  optimization  technique  of  taking  the  partial 

derivatives  of  Equation  (4.2)  with  respect  to  C  and  C, 

o  i 

and  setting  each  result  equal  to  zero.  The  resulting 


so-called  "normal  equations"  are  then  solved  simultane¬ 
ously  (Neter  &  Wasserman,  1974) . 
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If  q  =  1  in  Equation  (4.1),  the  solution  technique 
for  the  generalized,  bivariate,  linear  location  problem 
is  not  as  straightforward.  The  following  result,  after 
the  substitution,  prohibits  the  use  of  the  classical 
optimization  technique: 


CAl-<Cc 


cl>  = 


-  C. 


"  Claill 


(4.3) 


One  of  the  first  solution  techniques  used  to  solve 
Equation  (4.3)  was  proposed  by  Chames,  Cooper,  and 
Ferguson  (1955) .  They  demonstrated  that  the  following 
linear  program  produced  an  optimal,  although  not  neces¬ 
sarily  an  unique,  solution  (Sposito &  Smith,  1976). 

Let  dt  and  d^  be  the  vertical  deviations  "above"  and 
"below"  the  linear  facility  to  be  located  for  the  ith 
existing  facility.  Thus,  for  any  existing  facility  i 


C..a. ,  -  C  +  d.  -  d.  =  a. 

1  il  o  i  l  i2 


The  simplex  formulation  of  the  problem  takes  on  the 
following  structure : 


Minimi  ze 


E  (W.d+  +  W.dT ) 
i=l  “  1  11 


(4.4) 


subject  to 


Clail 


+  C  +  dT 
o  i 


>_  0 ,  and 

are  unrestricted. 


i 2 


(4.5) 


With  modifications  for  unrestricted  variables,  any  simplex 
program  can  be  used  to  solve  this  linear  programming 
formulation. 

Linear  programming  can  be  used  to  solve  another 
member  of  this  family — specifically,  when  q  approaches 
infinity  (Chebyshev  criterion) .  When  q  approaches 
infinity,  the  problem  becomes  the  location  of  the  linear 
facility  such  that  the  maximum  vertical  distance  from  the 
existing  facility  to  the  linear  facility  is  minimized. 
Wagner  (1959)  demonstrates  the  formulation  which  follows. 
Again,  let  d  represent  the  deviations  from  the  existing 
facility  to  the  linear  location  line  to  be  located.  In 
this  case,  however,  let  d  be  common  to  each  of  the  follow¬ 
ing  inequalities: 


-d  <  C1ail  + 


Co  - 


ai2  i  +d 


V 


(4.6) 


Therefore,  the  linear  program  becomes 
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Minimize  d 
subject  to 


C,a.  ,  +  C  - 
1  ll  o 

d  i  ai2 

Vi 

C, a. .  +  C  + 
1  ll  o 

d  i  ai2 

V. 

1 

d  >  0 ,  and 

C,  and  C  are  unrestricted.  (4.7) 

1  o 

Again,  any  simplex  program  can  be  used  to  solve  the  pre¬ 
ceding  formulation  (i.e.,  with  slight  modification  for 
unrestricted  variables  depending  upon  the  program) . 

Forsythe  (1972)  extended  the  range  of  analytical 
solutions  for  this  family  for  values  of  1  <  q  <  2  by  the 
use  of  gradient  search  techniques,  and  Morris  and  Norback 
(1980)  extended  the  range  for  values  of  q  between  2  and 
infinity  by  the  use  of  convex  programming  and  Kuhn-Tucker 
conditions.  The  applicability  of  these  techniques  results 
from  the  property  which  follows . 

Property  4.1 

CLq(C  ,  C. )  is  convex  with  resptct  to  both  C  and  C. 

p  O  1  O  1 

when  the  limit  is  taken  as  p  approaches  1  from  either  the 
positive  or  the  negative  direction  and  q  1. 


6 


T 


Proof.  The  Weber  formulation  with  the  hyperbolic 
smoothing  constant  for  the  generalized,  bivariate,  linear 
location  problem  when  the  limit  is  taken  as  p  approaches 
1  from  the  negative  direction  is  as  follows: 


CLq_(CQ,  CL)  =  ,^iwinai2 -C0-C1ail)2  +  m]q/2. 


(4.8) 


The  first  partial  derivatives  of  CL^-(Cq,  C^)  with  respect 


to  C  and  C.  are  as  follows: 
o  1 


&CL?-(C  ,  C  )  n  2  (a/2)-l 

.  3C  *  E.V,ai2  -  Co  '  Vil>  +  ml 

O  1=1 


(ai2  ~  Co  “  Clail}  ' 


(4.9) 


3CLq-(Cp,  Cl) 


n 


3C. 


-qi^1MiHai2  -  Co  ‘  Clail>  + 


m] 


(q/2) -1 


*  (ai2  “  C~  “  C,a41)a41. 


Til  il‘ 


(4.10) 


The  second  partial  derivatives  of  CLq_(CQ,  C^)  with 
respect  to  CQ  and  are  as  follows: 


32CL?-(C  ,  C.)  n  . 

-  ae°2  -  ■J{JlV<ai2  -  Co  '  Clail>  +  ”> 


(q/2) -2 


*  l(ai2  "  Co  “  Clail)  (q  "  X)  +  ml}'  (4*11) 


I 


-  °o  -  w!  *  *> <q/2)-2 


*  t(ai2  “  Co  "  C1ail)  (q-D  +  (4.12) 


^  CL*  <C,  C,)  n  2  (cr/2) -2 

- -  q(  r AH.U  -  =0  -  V  >  +  m)(q/2) 

O  1  1=1 


*  f(ai2  "  Co  "  C1ail)^(q_l)  +  mJaii}*  (4-13) 

The  determinant  of  the  Hessian,  D,  for  the  family  of 
functions  where  the  limit  is  taken  as  p  approaches  1 
from  the  negative  direction.  Equation  (4.1),  is  as 
follows : 

32CL?_(C_,  C.)  32CL?_(C  ,  C.)  32CL?- (C  ,  C. )  , 

o  =  ( — y  f— -3 - )( — ^  -  ( — io? — ' 


-  q2({yw.[(ai2  -  Co  -  cl3il)2  +  m](q/2)"2 

*  t(ai2  "  Co  "  Clail)2(q  *  1}  +  mj} 


*  <J1Mit(ai2  -  Co  -  Clail»2  *  m]  *q/2’-2 

*  t(ai2  -  CQ  -  C1ail)2(q  -  1)  +  mja2^ 


-  -  Cc  -  Cllu)2  + 


*  t(ai2  "  Co  "  Clail)2(q  -  L>  +  ni]ail}2]. 


(4.14) 
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If  this  determinant  (4.14)  is  positive,  it  is  a  sufficient 
condition  for  convexity  since  both  Equations  (4.11)  and 
(4.12)  are  positive  when  q  1.  The  sign  of  Equation 
(4.14)  is  difficult  to  evaluate,  however,  because  of  the 
summations.  An  indirect  method  of  determining  the  con¬ 
vexity  properties  of  CL^_(CQ,  C1)  is  examination  of  the 

function  on  a  term-by-term  basis  (i.e.,  L?_.(C  ,  C, ) ) . 

1  l  o  l 

Given  that 

Ll-i(Co'  Cl)  =  I(ai2  "  C0"C1ail)2  +  mlq/2,  (4.15) 


the  relevant  partial  derivatives  with  respect  to  CQ  and 
are  as  follows: 


3Ll-i(Co'  ^*1^  2  (q/2)~^ 

-He  -  -S[<ai2  -  Co  -  Clail>  +  ml  q 


(ai2  "  Co  “  Clail) ' 


(4.16) 


9L?_.(C  ,  C.)  ,  f_/54  , 

Icf  - S‘<ai2  -  Co  -  clail>  +  »!  q/  ’ 


*  (ai2  -  Co  -  ClaU,ail' 


(4.17) 


32L?_. (C  ,  C.)  9  , 

■  l  f  -  ^<ail  -  co  -  Clail>  *  9  ( 

3Co 


*[(ai2  ”  Co  "  Clail)  (q  "  1)  +  m] ' 


(4.18) 
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32L?-i«V  Cl) 


3C, 


2 - *[(ai2  ”  Co  '  Ciail)2  +  m](<*/2)-2 

*  [(ai2  "  Co  "  Clail)2{q-1)  +m]ail2,  (4.19) 


32L?-i(CC  Cl-  -  ,t(ai2  -  C0  - 


3C  3C. 
o  1 


*  [(ai2  “  Co  '  C1ail)  (q-D  +mlaii*  (4.20) 


The  determinant  of  the  Hessian,  D,  for  Lq-^(Co,  c^)  is 
given  by  the  following  formula 


32L?_. (C  ,  C.)  32L?_ . (C  »  C.)  32L?_ . (C  ,  C. )  , 

D  =  ( - t-i— 2 - i-)  ( - ±-±—2 - — )  -  ( - 1  -1  ~° - M  , 


3C 


3C, 


3C  3C.. 
o  1 


-  12t  <ai2  -  Co  -  Clail'2  +  ■"l(q/2>'2 


*  I,ai2  •  Co  -  °lail,2,,5  -  11  +  -"Itai!2  -  au2). 


=  0. 


(4.21) 


Since  Equations  (4.18)  and  (4.19)  are  positive  when  q  >  1 
and  each  element  of  this  determinant  is  equal  to  zero, 
LjL^(CQ,  C^)  is  convex.  This  property  implies  that 
CLq_(C0,  Cx)  is  convex  since  it  is  the  sum  of  convex 
terms,  and  the  property  is  proven. 

For  0  <  q  <  1  Equations  (4.11)  and  (4.12)  are  posi¬ 
tive  indicating  CLjL(CQ,  C^)  is  concave.  Therefore, 
gradient  search  techniques  and  convex  programming  are  not 
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conducive  to  determining  optimality.  Barrodale  and 
Roberts  (1970)  proved  that  the  values  of  CQ  and  which 
optimize  Equation  (4.1)  when  0  <  q  <  1  must  pass  through 
at  least  two  points  (a.^,  akl'  ak2  wliere  i  i*  k)  . 

Thus,  a  simplistic  solution  is  the  determination  of  each 
possible  line  that  passes  through  two  data  points  and 
determination  of  the  value  of  C£q_(CQ,  C^)  for  each  line. 
(See  Barrodale  &  Roberts,  1970,  for  a  proof  of  the  above 
procedure.  A  variation  of  the  proof  is  presented  later 
in  this  chapter.) 


The  p  =  2  Family 

The  Weber  reformulation  of  the  p  -2  family  of  the 
generalized,  bivariate,  linear  location  problem  is  given 
by  the  following  formula  which  is  a  special  case  of  Equa¬ 


tion  (3.13) : 


“I'V  C1»  - 


i;1wil<»i2  -  co  -  cian)2  + m|q/2 

[i  +  c12]q/2 


(4.22) 


The  members  of  this  family  of  the  generalized,  bivariate, 
linear  location  problem  seek  values  of  CQ  and  that 
minimize  the  perpendicular  distance  from  the  existing 
facilities  to  the  linear  location  facility  taken  to  the 
q  power.  Although  the  publications  in  this  area  have  not 


I 


l 
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been  as  plentiful  as  those  in  the  p  *  1  family,  extensive 

research  has  been  done  for  certain  cases. 

The  q  =  2  model  is  the  most  common  application  of 

the  p  =  2  family.  More  commonly  known  as  factor  analysis, 

this  model  determines  values  for  C  and  C.  that  minimize 

o  1 

the  expression 


n 

I 

CL2<Co'  °1>  *  — 


WiIai2  *  Co  *  Clail! 


U  +  C 


(4.23) 


Classical  optimization  cam  be  used  to  solve  the  preceding 
expression  even  though  the  typical  method  employed  involves 
formulating  the  eigen- structure  (Green,  1978) .  In  the 
simple  bivariate  case,  the  following  procedure  is  adequate. 
1.  Determine  the  first-order  partial  derivatives: 

8CL2(Co-  Cl> 

3Co 


■2 »i <ai2  -  Co  -  Clail> 
l—i  2 

[1  +  ] 


(2.24) 


dCht (C  ,  C.) 


2  g-"1-  *  -2  {.^WiIai2  “  Co  -  Clail][1  +  Cl  ]a 


il 


+  ci  .£^12  -  co  -  ciaii]  >• 


(4.25) 
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2.  Set  Equations  (4.24)  and  (4.25)  equal  to  zero, 
and  solve  simultaneously: 


c  -  -A  i  (A2  +  4B2)1/2 
U1  “  2B 


(4.26) 


where 


n  n 

A  =  E  W. [  I  W.a  ‘ 

i=l  1  i=i  1 


.Vu2i  +  ‘ ^u)2 

1=1  1=1 


-  (  £  W  a.  )  , 
i=l  1 


n 


n 


n  n 


B  =  E  W.a.  E  W. a . -  -  E  W.  E  W  a  a.  ,  and 
i=l  ^  ^  i— 1  ^  i=l  ^  i— 1  ^ 


C  = 


.^Wi(ai2  "  clail] 
n 

i£i"i 


(4.27) 


The  preceding  formulation  will  result  in  two  solutions: 
maximizing  Equation  (4.23)  and  minimizing  Equation  (4.23). 

Morris  and  Norback  (1980)  extended  the  work  of 
Wesolowsky's  numerical  technique  into  an  analytic  solution 
for  Equation  (4.22)  when  0  <  q  <_  1.  Using  a  technique 
similar  to  that  of  Barrodale  and  Roberts  (1970)  for  the 
vertical  distance,  Morris  and  Norback  proved  that  a  line 
minimizing  Equation  (4.22)  must  also  pass  through  at  least 


two  existing  facilities.  Thus,  the  solution  technique 
for  0  <  q  £  1  and  p  =  2  is  identical  to  the  technique  used 
to  solve  the  limit  as  p  approaches  1  from  the  negative 
direction  family. 

Morris  and  Norback  pointed  out,  in  the  same  article, 
that  the  techniques  used  to  solve  the  vertical  distance 
problem  for  q  >  1  cannot  be  used  for  the  p  =  2  family 
since  Equation  (4.22)  becomes  nonconvex  and  nonconcave 
when  q  >  1. 

The  g  =  2  Family 

The  q  =  2  family  of  the  generalized,  bivariate, 
linear  location  problem  is  significant  for  two  important 
reasons:  (1)  linear  regression  (p  approaches  1  and 
q  =  2)  and  factor  analysis  (p  =  2,  q  =  2)  are  the  most 
widely  used  applications  of  linear  models  today,  and  (2) 
the  mean  which  is  the  basis  of  the  q  -  2  model  allows  for 
simplified  solution  techniques  that  provide  analytical 
results  regardless  of  the  p  value.  Since  parametric 
statistics  are  based  upon  the  mean,  the  q  =  2  model  has  a 
very  strong  theoretical  attractiveness  even  though  it  is 
overly  sensitive  to  outliers.  The  simplified  solution 
techniques  for  the  q  =  2  model  result  from  the  property 
which  follows. 
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Property  4 . 2 

Given  the  function 


=  "!“•«  -  Co  -  Vu>2  *  -"19/2 

CL?<Co'  Cj)  =  1-1 


[1  +  ((c12)1/2]p/I>'i]'J(p-l,/P  ’ 


(4.28) 


*  *  2 

the  values  of  C  and  C.  that  minimize  CL  (C  ,  C.. )  result 
o  1  p  o  1 

in  a  line  that  passes  through  the  point  (a^,  a2)  where 


AWiail  _  ="1*12 

‘1  =  “n -  and  a2  =  ^ - 


I  W 
i=l  x 


z  w 
i=l  1 


Proof .  If  q  =  2,  then  Equation  (4.28)  reduces  to: 


n 


clUc.  C.)  = 
p  o  1 


l  W. (a . _  -  C  -  C.a.  ) 
1  12  O  1  U 


i=l 


u  +  [(c2)l/2jp/(p-l)]2(p-l)/p* 


(4.29) 


*  _  ★  _ 

Let  a^  =  a^  -  a^  and  a^2  =  a^2  -  a2  be  the  mean  cor¬ 
rected  values  for  the  existing  facilities  (a^,  a^2)  * 

*  * 

Substituting  the  values  a^  and  a^2  into  Equation  (4.29) 
and  taking  the  first  partial  derivative  with  respect  to 


CQ,  the  following  result  is  obtained: 
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3L^(Co,  CL)  "2i=1Wi(ail  C°  Clail) 

3C0  =  “  ^(c2)l/2]p/(p-l)]2(p-l)/p-  <' 

Setting  Equation  (4.30)  equal  to  zero  and  substituting 
★  _  ★  _ 

a^  =  a^  -  a^  and  a ^  “  a2  results  in  the  following 
equation: 


i^lWi(ai2  a2  “  Co  Cl(ail  “  al))  "  °' 


n 


i-1 


w.  a .  _ 
i  x2 


n  ___ 

-  £  W.a-  - 

i-i 1 2 


n  n 

I  W.C  -  C,  £  W.a. , 

i=i  1  °  k-i  1  11 


and  therefore. 


-  ail  -  °- 


Cc  •  0.  (4.31) 


By  mean  correcting  the  data,  the  coordinates  are  trans¬ 
lated  such  that  the  point  (iu  ,  a2>  is  located  at  (0,  0) . 
Since  CQ=  0  when  the  data  is  mean  corrected  regardless  of 
the  p  value,  the  property  is  proven. 

With  the  knowledge  that  =  0  when  q  =  2  and  the 
data  is  mean  corrected,  Equation  (4.29)  can  be  rewritten 
as  follows: 


I 
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clUc . ) 
p  1 


n  * 

.*  "ita« 
1  =  1 


-  c^)2 


(4.32) 


Taking  the  first  derivative  with  respect  to  C-L  and  setting 
it  equal  to  zero  produces  the  following  results: 


dCIT(C, ) 

_E _ i-  - 


n 

-2  2  w. 

i=l  z 


dC, 


U  7 

*  ([1  +  UC^|1/2(p/p-1,l <a*2 
+  c1nc2i1/2)(2-p)/,p-1)  C«u 


claU,all 


claU,K 


(4.33) 


[(c2)l/2,(p-2)/(p-l)  j  „.(a*2  - 
x  i=l 


*  * 


c2  2  W. (a,„  -  Ciail)au 


1i=l  1  i2 


n  *  *  2 

+  Cll=1MiUi2  '  Clail’ 


-  0, 


(4.34) 


n 


2  W. [a.2 

[(c2, 1/2,  1/<P-D  ({i=i - - 


*  3  * 

Clail  ai2 


E  Witai2 
i=l 


Claillail 


]2}l/2 


1.  (4.35) 
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Let  p  be  equal  to  a  rational  number  u/v.  Then  Equation 
(4.35)  can  be  reduced  to: 


2v, 


n 

* 

*  * 

=  Wi 

Ui2 

'  Clail>ai2 

i=l 

n 

* 

★  ★ 

z  w. 
i=l  1 

(ai2 

Clail)ail 

=  0. 


(4.36) 


Given  the  assumption  that  p  is  a  rational  number 
(i.e.,  not  an  overly  restrictive  assumption).  Equation 
(4.36)  can  be  solved  for  any  value  of  p.  If  the  resulting 
polynomial  is  greater  than  order  4,  a  numerical  technique 
such  as  the  Birge-Vieta  iterative  method  (McCalla,  1967)  is 
required  to  determine  all  possible  solutions  for  C^.  The 
value  of  that  minimizes  Equation  (4.32)  can  then  be 
found  by  straight  substitution.  Theoretically,  the  pro¬ 
cedure  outlined  in  this  section  is  applicable  for  any  values 
of  p  (i.e.,  as  long  as  p  can  be  considered  a  rational 
number) .  Certain  values  of  p,  however,  result  in  polynomials 
of  excessively  large  order  (i.e.,  p  =  1.51  produces  a  301th 
order  polynomial,  while  p  =  1.50  produces  only  a  6th  order 
polynomial) .  In  most  cases ,  use  of  the  heuristic  solution 
presented  later  is  easier. 

The  q  =  1  Family 

The  q  *  1  family  of  the  generalized,  bivariate, 
linear  location  problem  describes  the  median  line  that 


2 


minimizes  the  absolute  distance  in  the  p-norm  direction. 
The  q  =  1  family  has  an  advantage  in  that  the  located 
linear  facility  is  less  susceptible  to  outliers.  For 
q  =  1  the  objective  function  to  be  minimized  is  given  by 
the  following  equation: 


CLp(Co' 


Cl> 


I  W.|a.  -  C 

i=l  1 


-  C.a.J 
o  1  ll1 


[i  +  t(c^)1/2]p/(P"i)]  (P“D/P 


(4.37) 


The  following  property  implies  a  simplified  solution 
technique  for  the  q  =  1  family. 

Property  4.3 

*  * 

A  solution  (CQ  and  C-^)  to  Equation  (4.37)  exists 

it  fc 

which  satisfies  a ^  =  CQ  +  ciaii  for  at  least  two  existing 

facilities  (i  and  k  where  i  /  k). 

The  following  proof  is  a  generalization  of  the  method 

used  by  Morris  and  Norback  (1980)  for  the  p  *  2  case. 

Proof.  Since  a  solution  of  the  form  x_.  ^  =  C_  +  C,x.., 
-  12  o  1  ll 

exists,  the  axes  can  be  translated  without  a  loss  of 
generality  so  that  CQ  =  0.  Then  CJlp(C1)  can  be  written 
as  follows: 


c*p(ci>  ‘  tl 


(4.38) 


a 
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The  existing  facilities  i  can  be  renumbered  so  that  the 
following  region  is  specified: 

Rr  =  ai2  “  Clail  —  i=l»  ...,r; 

ai2  “  Clail  -  °'  1  *  r  +  1'  *•*'  n}*  (4*39) 


The  union  of  regions  R^.,  r  =  0,  . ..,  n  covers  the  real 

number  line.  When  r  =  0  (r  =  n)  no  values  of  i  exist 

such  that  ai2  -  C^a^  :>  0  (ai2  -  c1ail  <  °)  •  For  e  Rr 

the  absolute  value  sign  can  be  removed  from  0.^(0.)  and 

P  1 

the  following  equation  can  be  written: 


CJlptC^  -  [1  +  [(C2)1/2]p/(P“1,]~(p“1)/p 


*{i=lWi(ai2  "  Vil* 


Z  W.(a  -Ca  )}, 
i=r+l  1  1  1X 


-  [1+  [  (Cj)1/2]P/(P"1,]“(P”1)/P(A-C1B) 

(4.40) 


where 


A  =  Z  W.a.0  -  Z  W.a.,, 


t-i  1  12  i 


i=r+l 


i  121 


£  ^iail  ~  ^  ^iail' 

i=l  1  11  i*r+l  1  11 


B 
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The  term  (A  -  C^B)  in  Equation  (4.40)  must  be  greater  than 

zero  unless  the  existing  facilities  are  collinear,  in 

which  case  the  proof  is  trivial. 

All  that  remains  to  be  shown  is  that  cannot 

P  1 

be  minimized  on  the  interior  of  Rr>  The  first  derivative 
of  C£p(C1)  is  as  follows: 


dCJr(C..) 
_ P _ 1 


dC, 


■  -Cl  +  [(c2)1/2]_[(2p~1)/pl] 

*  (C1f(C2)1/2]  (2_p)/(p''1)  A+  B)  . 


(4.41) 


Setting  Equation  (4.41)  equal  to  zero  results  in  the 
following  equations: 

(4.42a) 

(4.42b) 


or 


[ (C2)1/2] (2-p)/(P"1)C1  =  -[B/A] 
=  ± [ (B/A) 2] (p“1)/2. 


Taking  the  second  derivative  of  Equation  (4.41)  produces 
the  following  derivative: 


d2ciii  (cn ) 
_ E _ L. 


[(c2)l/2] [p/(p-l)]-2 

(FT) 

U  +  nc2)l/2]p/(p-l)]-[(3p-l)/p] 
[2(p-l) [ (C2)1/2]p/ (p-1)A  +  (2p-l) 


*  C.B  -  A] 


3 
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The  sign  of  d2^^)  /dC2  is  determined  by  the  expression 

[2(p-l)[(C2)1/2]p/(p‘1)A+  (2p-l)C1B-A].  When  the 

values  of  C,  from  Equation  (4.42b)  are  substituted  into 

2  1  2 

the  preceding  expression,  d  C£p(C^)/dC^  has  the  same 
sign  as 

2 (p  -  1) [ (B/A) 2]P/2A  ±  (2P  -  1) I (B/A) 2)  (P_1) /2B  -  A 


s  Cr(B/A)21 (P"2)/2B2  +_A2} 
A 


(4.44) 


Thus,  d2CAp(C1) /dC2  has  the  same  sign  as  -A.  A  -  C^B  >  0 
j_g  violated,  however,  when  takes  on  the  sign  from 
Equation  (4.42a)  unless  A  >  0.  Therefore,  any  extremum 
on  the  interior  of  Rr  must  be  maximum  when  p  ^  1.  Since 
CAp(Ci)  is  bounded  below  on  Rr,  CAp(C^)  is  minimized  on 

the  boundary  of  R^.  where  at  least  one  term  of  the  form 

* 

Wj:  aj2  '  °laji  =  °- 

Substituting  C1  in  Cip(C0,  Cj) ,  the  condition  func¬ 
tion  can  be  written  as 


c*p<cclcl>  - 


x4+iWi'a« '  c°! 


(4.45) 
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where 

• _ ^ _ 

Wi  [l  +  j(C  )l/2jP/(p-l) j (P-1)/P# 

ai2  ~  ai2  "  Clail'  and 

the  jth  term  of  C£^(C  ,  C, )  which  equals 

p  o  l 

zero  is  omitted. 

is  convex  from  Property  3.4  and  piecewise 
linear  with  points  of  discontinuity  in  the  first  deriva- 

f 

tive  occurring  at  a^/  i  =  1/  • • • *  n*  i  £  j •  A  mini¬ 
mizing  value  of  CQ  exists  for  which  at  least  one  term  of 
the  form  W.  |  a.  2  “  C  |  “  0  and  j  t  k.  This  condition 

means  that  at  least  the  jth  and  kth  terms  of  C£1(C  ,  C. ) 

p  o  1 

equal  zero,  and  the  property  is  proven  for  all  p  except 

P  =  1. 

The  results  of  the  previous  property  extend  the 
simplistic  technique  developed  by  Morris  and  Norback 
(1980)  for  p  =  1  and  q  =  1  to  the  general  case  of  any 
value  of  p.  To  solve  any  q  =  1  generalized,  bivariate, 
linear  location  problem,  determine  the  line  that  passes 
through  two  existing  facilities  and  minimizes  C£p(CQ,  C^). 

The  p  «  <*>  Family 

The  p  =  *  family  of  the  generalized,  bivariate, 
linear  location  problem  is  the  special  case  where  the 
error  is  divided  equally  between  the  two  variables 
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and  When  p  =  00  (i.e.,  as  with  all  cases  when 

p  5*  1)  ,  the  concept  of  independent  and  dependent  variables 
no  longer  holds.  Each  variable  in  the  p  *  00  model  must 
be  responsible  for  an  equal  share  of  the  error  introduced 
by  the  linear  model.  When  p  =  00  the  generalized,  bivari¬ 
ate,  linear  location  model  becomes 


n 


^1*12  -  Co  -  Vill 

U  +  |c1|]<3 


(4.46) 


The  distance  being  considered  in  this  family  is  the 
45-degree  distance  from  the  given  facility  to  the  linear 
facility  to  be  located.  This  result  is  derived  from 
Property  2.7  where  the  statement  is  made  that  the  slope 
of  the  line  between  the  given  facility  i  and  the  linear 
facility  to  be  located  that  minimizes  the  p-norm  distance 
is  given  by  1/[(C2)^2]^^  When  p  -  °°  the  slope  is 

equal  to  1.  The  Weber  formulation  of  the  p  =  »  family 
with  the  hyperbolic  smoothing  constant  is  as  follows: 


cl?  » 


Z  W,[(ai2  -  C0  -  Cja^)2  +  m]q/2 


[1  +  (c2)1/2]q 


(4.47) 


The  special  case  of  p  =  «  and  q  *  2  is  interest¬ 
ing  because  of  its  close  relationship  to  regression.  If 


If  Equation  (4.35)  is  rearranged,  the  following  equation 
is  obtained: 


If  p  approaches  1  from  the  negative  direction  as  in  the 
linear  regression  of  x^  on  x.^,  the  limit  of  the  term  on 
the  right  side  of  the  equality  approaches  infinity.  The 
necessary  and  sufficient  condition  for  the  equality  to 
hold  is  for  the  denominator  of  the  term  on  the  left  side 
of  the  equality  to  approach  zero.  The  denominator  of  the 
term  on  the  left  side  of  the  equality  set  equal  to  zero  is 
the  normal  equation  for  linear  regression  with  mean- 
corrected  data,  and  the  term  to  the  right  of  the  equality 
is  the  slope  of  the  line  connecting  the  existing  facility 
i  to  the  point  on  the  linear  facility  that  minimizes  the 
p  distance. 

A  similar  relationship  holds  when  p  approaches  1  from 
the  positive  direction.  The  limit  o'  the  term  to  the 
right  of  the  equality  is  zero  as  the  limit  is  taken  as 
p  approaches  1  from  the  positive  direction.  The  left  side 
of  the  equality  can  be  equal  to  zero  only  when 
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n 

E  W.  [a 
i=l  1 


* 

i2 


cia‘n 


=  0. 


As  shown  by  Property  3.6,  the  preceding  expression  is  the 
normal  mean-corrected  equation  for  the  linear  regression 
of  on  X2» 

If  the  limit  of  the  right  side  of  the  equality  in 
Equation  (4.35)  is  taken  as  p  approaches  infinity,  the 
following  relationship  results: 


[  E  W.  (a 
i=l  1 


* 

i2 


Clalllai2i 


2 


n 

=  [  E  W. [a 
i=l  1 


* 

i2 


Clai2 


*  *  2 
Clail]ail]  ' 


(4.49) 


The  preceding  equation  results  directly  from  the  fact 
that  the  limit  of  1/  [  (C^) 1^2]  is  1  as  p  approaches 

infinity.  Equation  (4.49)  can  be  solved  directly  by 
expanding  terms  and  using  the  quadratic  formula.  Thus, 
Equation  (4.49)  can  be  considered  the  normal  equation  for 
p  =  oo ,  q  =  2 . 

Note  that  values  of  p  approaching  1  from  the  positive 
or  negative  direction  and  of  p  approaching  infinity  are 
the  only  cases  of  the  generalized,  bivariate,  linear 
location  problem  where  1/ [  (C^) 1^2]  is  a  constant. 

For  all  other  values  of  p,  the  line  connecting  the 
existing  facility  i  to  the  point  on  the  linear  facility 
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that  minimizes  the  p  distance  is  a  function  of  the  linear 
relationship  (C^) . 

The  0  <  q  <  1  Family 

The  0  <  q  <  1  family  of  the  generalized,  bivariate, 
linear  location  problem  deals  with  the  general  case  where 
the  existing  facilities  closest  to  the  linear  facility  to 
be  located  have  the  greatest  influence  in  the  decision 
process.  The  problem  of  determining  the  location  of  the 
linear  facility  poses  several  difficulties  in  determining 
optimal  solutions  because  of  the  nonconvexity,  the  non¬ 
concavity,  and  the  nondifferentiability  of  the  objective 
function.  Certain  simplifications  are  available,  however, 
that  ease  the  computation  difficulty. 

As  stated  in  Property  3.4,  Lp^(CQ,  C^)  is  strictly 
concave  with  respect  to  CQ  for  0  <  q  <  1  since  the  hyper¬ 
bolic  smoothing  constant  m  is  greater  than  zero. 

C£jj(C  ,  C.)  ,  however,  contains  discontinuities  in  the 

p  O  X 

first  derivatives.  These  discontinuities  allow  for  the 
proof  which  follows . 

Property  4.4 

Given  the  following  formulation  of  C£q(C  ,  C. )  : 

p  o  l 

CA?(C  ,  c.) 
p  o  1 


n 


q 


-  Co  -  Clail 


(4.50) 
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an  optimal  minimizing  solution  exists  for  CQ  and  called 

*  *  *  * 

C  and  C,  which  satisfies  a.~  =  C  +  C.a..  for  at  least 
o  1  i2  o  1  il 

one  existing  facility  i  when  0  <  q  <  1. 

Proof.  For  a  proof  by  contradiction,  first  assume  a 

*  * 

solution  exists  for  C  ,  C.  called  c  and  C,  which  mini- 

o  1  o  1 

mizes  Equation  (4.50)  and  which  does  not  satisfy 

4r  k 

ai2  =  CQ  +  C1ail  for  any  existing  facility  i.  Then  the 
existing  facilities  can  be  partitioned  into  two  subsets 

•k  k 

depending  upon  whether  a^2  “  CQ  “  ciaii  *s  <Jreater  than 
or  less  than  zero.  C&p(C*,  C*)  can  then  be  rewritten 
within  the  absolute  value  signs  as  follows: 


Li» i,ai2  -  Co  -  Vil)q 


c  )  —  ^  - 

C  p(Co'  Cl'  ■  fl  +  [(c*2)l/2]p/(p-l)]q(p-l)/p 


n 


i=LlWi(C°  +  Clail  ~  ai2} 

[i  +  [  (c*2)1^2]p/,^p-1)  ]q^p_1)/,p 


(4.51) 


where 

0  <  q  <  1. 

* 

Given  the  value  of  the  slope  that  minimizes  Equation 

(4.50),  the  conditional  first  and  second  derivatives  for 

CAq(C*|ct)  are  as  follows: 
p  o '  1 
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dC£?(C*|  C.  ) 


=  -q  [  1  +  nc*-2)l/2]p/(p-l)]-q(p-l)/p 
dC  1 


*  *  q-1 


*  {/=1Wi(ai2  "  Co  -  Clail) 


n 


*  *  cr—  1  ■, 

E  W.  (C  +  C.a  -  a .  _ ) 4  , 

i-k+1  1  °  1  11  12 


(4.52) 


dCo 


Z  W. (C*  +  C*a..  -  a.-)q~2,.  (4.53) 

i=k+l  1  °  1  11  1Z  ' 


The  second  derivative  Equation  (4.53)  is  strictly  less 

★ 

than  zero  when  q  <  1.  Thus,  CQ  corresponds  to  a  relative 
maximum  rather  than  to  a  relative  minimum.  This  condition 
requires  that  the  minimizing  linear  facility  pass  through 
at  least  one  existing  facility,  and  the  proof  is  complete. 

The  alternate  proof  mentioned  earlier  for  the 
Barrodale  and  Roberts  (1970)  theorem  that  the  linear 
facility  minimizes  the  vertical  distance  where  0  <  q  <  1 
passes  through  at  least  two  existing  facilities  is 
presented  here. 
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Property  4.5 


Given  the  following  formulation  of  c£jL(CQ,  C^) : 


C*l“(Co'  Cl)  =  "  Co  '  ClaiJq 


(4.54) 


where 

0  <  q  £  1, 

*  * 

an  optimal  solution  exists  for  CQ  and  called  CQ  and 
that  minimizes  C^)  and  satisfies  ai2  -  CQ  +  C^a^ 

for  at  least  two  existing  facilities  (  i  —  j  and  i  =  k  and 
j  /  k)  when  0  <  q  <  1. 

Proof .  From  the  results  of  Properties  4.3  and  4.4 
the  statement  can  be  made  that  the  linear  facility  that 
minimizes  CZ^-(Co,  C1)  will  pass  through  at  least  one 
existing  facility.  Showing  that  the  minimizing  linear 
facility  must  pass  through  a  second  existing  facility 
completes  the  proof.  Using  a  method  similar  to  the  one 
used  for  Property  4.4,  translate  the  coordinates  so  that 


the  one  existing  facility  (a^,  that  is  on  the  mini¬ 

mizing  linear  facility  is  at  the  origin.  C£^L(C0,  C^) 


can  now  be  rewritten  as: 


(4.55) 


where 


i  =  1,  ...»  k-1  are  the  existing  facilities  that 
satisfy  ai2  “  cian  > 

i  =  k  +  1,  . ..,  n  are  the  existing  facilities 
that  satisfy  C^a^  -  ai2  >  0,  and 

Wk(ak2  "  Clakl>q  =  °* 

* 

If  the  assumption  is  made  that  the  value  of  C1  called 
for  the  linear  facility  that  minimizes  Equation  (4.55) 
does  not  pass  through  a  second  existing  facility,  then 
C£q_(C^)  is  both  continuous  and  dif f erentiable .  Straight 
forward  differentiation  of  Equation  (4.55)  results  in  the 
following  equation: 


dC*q  (C.)  k-1 

- —  —  ■  —  =  -q  [  Z  W 


i(ai2  -  Clail) 


i=k+l 


Wi(Clail 


-  a  )q_1l 
ai2' 


(4.56) 


d2CHq-(C. )  k-1  _  - 

- i-T-i-  -  q(q  -  1)  [  I  Wi(a.2  - 


+  Z  W.  CC, a .  1  -  a.,)4* 
i=k+l  1  11  “ 


(4.57) 


Since  the  second  derivative  is  strictly  negative  for 

* 

0  <  q  <  1,  it  follows  that  any  must  be  a  relative 


maximum  rather  than  a  relative  minimum  as  assumed.  This 
contradiction  requires  that  the  relative  minimum  pass 
through  a  second  existing  facility,  and  the  property  is 
proven.  (The  proofs  of  Properties  4.4  and  4.5  are  based 
on  a  method  from  Morris  &  Norback  (1980)  suggested  by 
an  anonymous  referee . ) 

Property  4.5  can  be  easily  extended  to  the  limit  as 
p  approaches  1  from  the  positive  direction  by  simply 
reversing  the  variables  of  the  original  problem.  Thus, 
the  horizontal  distance  of  the  original  problem  becomes 
the  vertical  distance  of  the  modified  problem.  This 
process  is  similar  to  rotating  the  axis  of  the  original 
problem  by  90  degrees.  Property  4.5  was  extended  to  the 
p  =  2  case  by  Morris  and  Norback  (1980)  when  they  rotated 
the  original  axis  such  that  the  optimal  linear  facility 
was  horizontal.  In  the  p  =  2  case  the  vertical  distance 
is  the  distance  to  be  minimized  since  the  Pythagorean 
theorem  holds.  For  p  5*  2  or  p  f  1  the  rotation  method  is 
not  a  valid  solution  technique  because  of  the  triangular 
inequality  (Morris,  1981).  However,  the  following 
property  generalized  from  Morris  and  Norback  (1980)  eases 
the  computational  problems  in  determining  an  optimal  solu¬ 
tion  for  p  j*  2  or  p  t  1. 
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where 

C*  ,  I C? I  are  the  intercept  and  the  absolute 
ov'  *1 ' v 

value  of  the  optimal  slope  of  the  linear 
facility  when  p  approaches  1  from  the 
negative  direction, 

C*  ,  1 C? I  are  the  intercept  and  the  absolute 
op'  1  l'p 

value  of  the  optimal  slope  of  the  linear 

facility  for  any  p,  and 

C*.  ,  |C*|.  are  the  intercept  and  the  absolute 
oh  l  h 

value  of  the  optimal  slope  of  the  linear 
facility  when  p  approaches  1  from  the 
positive  direction. 

proof.  From  Equations  (3.15),  (3.16),  and  (4.58) 
the  following  equations  can  be  stated: 
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ciq<c,  c.)  »  ti  +  jc  |p/(p"1)i”q(p"1)/p 

P  O  J.  1 


*  Ctf.(C0,  Cj), 


(4.59] 


Ctg(C0,  Cl>  =  Ic.lu  +  |cl|P/<P-ll)-<5‘P-1,/P 


*  =‘i+(Co'  cl>- 


(4.601 


NOW  Cl’-(C*V,  c*v)  <  Cif_(C*p,  C*p>  since  (COT,  C'„) 
minimizes  CJljL(C0,  C^).  So  tl+|CjJp^*p  )  q^p  ^  ^p 

*  «?-<<„.  c*iv>  i » ♦  ic*li?/<p-urq,p-ll/E>ci?.(c;p,  c*p> . 
From  Equation  (4.45)  and  the  previous  expression, 

«J(C^.  C*lp)  <  Ci|<<v.  C*y)  -  <1  ♦  |C*|P/,p-1'rq,P‘1)/P 

*  CAq_(C*v,  C*v) .  This  means  [1  +  |C^|p/ (p_1) ]  q(p  1)/p 

*  11  ♦ 

so  that  |C*|P/(P”1)  1  |C*|P/(P_1).  Similarly, 

|c^hti+  lc*1lg/{p“1,l-q(p'1)/pc^(c^  <1>  i  l<lh 

*  11  +  |Cilg/(p"1)]'q(p-1)/pCiq+(C*p,  c*p)  and 

Icjlpti  +  |cJ|J/lp,"1)i"q<p’1)/pc*;+(cop,  clp>  <  |<lh 

*  [1  +  |C^|^/(p"1)]”q(p“1)/pCiJ+(C*h,  C*h) .  Therefore, 

icjipd  ♦  ic*ip/(p-1>rq(p"1)/p  <  i<ih 

*  [1  +  lci |p/,^p"1h-q(p"’1^p  which  resolves 

jc*  ^/(p-i)  <  jc*|p/(p-l)  completes  the  proof. 
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A  general  Heuristic  solution  for  the  0  <  q  <  1 

family  is  to  individually  select  each  existing  facility 

as  the  optimal  CQ  (i.e.,  from  Property  4.4)  and  to  vary 
*  * 

between  C^v  and  for  the  slope  that  minimizes 
Equation  (4.58).  The  optimal  solution  will  be  the  exist¬ 
ing  facility  that  results  in  the  minimum  value  for  Equa¬ 
tion  (4.58) . 


The  q  >  1  Family 

When  q  >  1  and  the  distance  to  be  minimized  is 

neither  the  vertical  nor  the  horizontal  distance,  both 

C1^(C  .  C.. )  and  CLq(C  ,  C. )  remain  nonconvex  and  noncon- 
p  o  1  p  o  1 

cave.  The  only  solution  technique  available  under  such 
circumstances  is  a  two-dimensional  Heuristic  search. 

The  search,  however,  is  simplified  as  a  result  of  Property 

*  *  * 

4.6  since  the  range  for  must  be  between  Clv  and  Clh- 

The  suggested  procedure  is  to  step  through  values  of 

and  to  determine  at  each  step  the  optimal  value  of  CQ 

given  the  present  value  of  CQ  can  be  determined  by 

any  gradient  search  technique  since  Ci.p(CQ,  C^)  is  strictly 

convex.  The  optimal  values  for  CQ  and  can  then  be 

determined  by  finding  the  minimum  value  of  Ci^(CQ,  C^) . 

Sample  Problem 

To  compare  the  results  of  various  combinations  of  p 
and  q  values  for  the  generalized,  bivariate,  linear 
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location  problem,  a  sample  problem  selected  from  Neter 
and  Wasserman  (1974,  p.  94)  is  presented.  The  optimal 
values  of  CQ  and  determined  from  7  values  of  q  ranging 
from  .1  to  infinity  and  from  6  values  of  p  ranging  from 
the  vertical  distance  to  the  horizontal  distance  are  cal¬ 
culated.  The  data  points  used  are  as  follows: 


i: 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

ail: 

7 

4 

5 

1 

5 

4 

7 

2 

4 

2 

8 

5 

2 

5 

7 

1 

4 

5 

ai2S 

97 

57 

78 

10 

75 

€2 

101 

27 

53 

33 

118 

65 

25 

71 

105 

17 

49 

68 

W.: 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

The  fact  that  the  weights  for  each  existing  facility  are 
equal  to  one  does  not  restrict  or  limit  the  example.  The 
results  of  the  calculations  are  presented  in  Table  4.1. 
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CHAPTER  V 


IMPLICATIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  RESEARCH 

The  derivation  of  the  generalized,  bivariate,  linear 
location  problem  consists  of  two  subproblems.  The  first 
subproblem  involves  the  determination  of  the  point  on  a 
linear  facility  that  minimizes  the  p-norm  distance  to  an 
individual  existing  facility.  The  second  subproblem  con¬ 
sists  of  determining  the  optimal  linear  facility  that 
minimizes  the  sum  of  the  q  multiple  of  the  p-norm  distance 
from  all  the  existing  facilities  to  the  point  on  the  linear 
facility  determined  by  the  previous  step.  The  simultane¬ 
ous  solution  of  the  two  subproblems  results  in  the  fol¬ 
lowing  formulations 


C£;J(C  ,  C-) 
p  o  l 


n 


Z  W.  |  a 
i*l__ 


i2 


[l  +  |c1|P/(P"1)]q(p‘:L)/p* 


The  lack  of  convexity  of  the  resulting  formulation 
for  the  general  case  prohibits  a  universal  solution  tech¬ 
nique  of  the  generalized,  bivariate,  linear  location 
problem  for  all  combinations  of  possible  values  for 


p  and  q.  For  certain  combinations  of  p's  and  q's  an  exact 
solution  of  generalized,  bivariate,  linear  location  prob¬ 
lem  cam  be  determined.  For  the  remaining  combinations 
the  only  meams  of  determining  a  solution  is  through  a 
heuristic  approximation  procedure.  In  Table  5.1  the  solu¬ 
tion  techniques  available  for  each  combination  of  p's 
and  q's  aire  outlined. 

Implications 

The  generalized,  bivariate,  linear  location  problem 
expands  the  use  of  linear  location  theory  to  include  the 
concept  of  a  p-norm  metric  system  of  measurement.  In 
previous  linear  location  problems  only  the  vertical 
distance  (the  limit  as  p  approaches  1  from  the  positive 
direction)  and  the  Euclidean  distance  (p  =  2)  were  con¬ 
sidered.  The  generalized,  bivariate,  linear  location 
problem  allows  the  analyst  to  place  the  linear  facility 
while  considering  the  importance  of  each  coordinate  axis. 
If  the  coordinates  are  measured  in  different  units  pro¬ 
hibiting  conversion,  for  example,  minimizing  the  Euclidean 
distance  may  not  be  consistent  with  reality  since  the 
p  *  2  distance  assumes  an  inverse  importance  of  the  axis 
depending  upon  the  slope  of  the  located  facility.  By 
proper  selection  of  the  p-norm  distance  to  be  minimized. 


Table  5.1 


Values 


0  <  q  <  1 

q  -  i 

1  <  q  <  2 

q  -  2 

q  >  2 

2  <  q  <  00 

q 


Where 


Solution  Techniques  for  the  Generalized 
Bivariate,  Linear  Locat.'on  Problem 


Values  of  p 


p  =  1* 

msti 

MSS 

mUmm 

warn 

D,  F,  G 

F,  G 

D,  F,  G 

F,  G 

A,  D,  G 

D,  G 

D,  G 

D,  G 

C,  G 

G 

G 

G 

A /  B,  C,  E,  G 

E,  G 

B,  E,  G 

E,  G 

C,  G 

G 

G 

G 

C,  G 

G 

G 

G 

A,  C,  G 

G 

G 

G 

B,  E,  G 


Linear  Programming  Formulation 
Classical  Optimization 

Gradient  Search,  Convex  Programming,  and 

Kuhn-Tucker  Conditions 

Two  Point  Solution 

Birge-Vieta  Iterative  Method 

One  Point  Approximation  Solution 

General  Heuristic  Technique  (Approximation) 


*The  limit  as  p  approaches  1  from  either  the  positive 
or  from  the  negative  direction. 
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the  analyst  may  select  the  weighting  of  the  coordinate 
system  that  best  defines  the  importance  between  the  two 
axes . 

The  generalized,  bivariate,  linear  location  problem 
can  be  extended,  just  as  in  Euclidean  linear  location 
theory  and  regression,  to  cases  other  than  the  minimiza¬ 
tion  of  the  absolute  or  the  square  of  the  distance  being 
considered.  The  value  of  q,  which  is  the  multiple  of  the 
absolute  distance  between  the  existing  facility  and  the 
linear  facility  to  be  located,  can  be  varied  between 
greater  but  not  equal  to  zero  and  infinity.  For  the  most 
part,  the  solution  techniques  available  for  linear  regres¬ 
sion  extend  to  the  generalized,  bivariate,  linear  loca¬ 
tion  problem  when  q  is  less  than  or  equal  to  one.  In 
general,  however,  for  q  greater  than  one  the  generalized, 
bivariate,  linear  location  problem  must  depend  upon  a 
heuristic  search  procedure  that  does  not  guarantee  an 
optimal  solution. 

Recommendations  for  Further  Research 

Three  areas  should  be  considered  for  future  research 
concerning  the  generalized,  bivariate,  linear  location 
problem. 

1.  The  first  area  concerns  the  problem  of  deter¬ 
mining  the  optimal  combination  of  p  and  q  values  that 


should  be  used  in  any  one  problem.  This  problem  may  be 
of  a  subjective  nature  in  most  situations.  If  the  models 
presented  in  this  research  are  used  in  a  forecasting 
situation,  however,  the  problem  becomes  quite  objective 
since  the  desire  is  to  obtain  the  best  fit. 

2 .  A  second  and  obvious  extension  of  the  bivariate 
problem  is  the  multivariate  case.  A  more  obscure  exten¬ 
sion  is  the  multivariate  problem  where  the  p  values  for 
each  variable  are  not  equal.  The  problem  can  be  formu¬ 
lated  as  follows: 


*2' 
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3.  The  third  area  concerns  the  statistical  applica¬ 


tions  of  the  generalized,  bivariate,  linear  location 
problem.  Extensive  work  must  be  done  to  provide  sampling 
distributions  for  the  various  parameters  of  the  model. 
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